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Introduction 



The Chow form of an equidimensional quasi-projective variety is one of the basic objects of algebraic 
geometry and plays a central role in elimination theory, both from the theoretical and practical points 
of view. 

Let V C P" be an equidimensional quasi-projective variety of dimension r and degree D defined over 
Q. Its Chow form Ty is a polynomial with rational coefficients — unique up to a scalar factor — 
which characterizes the set of over-determined linear systems over the projective closure V . More 
precisely, let Uq, . . . ,Ur denote r+1 groups of n -f 1 variables each , and setL^ := Uio xo + ■ ■ ■ + UinXn 
for the linear form associated to the group Ui ioi < i < r. Then G Q[f^o, ■ ■ ■ , Ur] is the unique 
— up to a scalar factor — squarefree polynomial such that 

(mQi • ■ • I "r) = V n {Lq{uq, x) — 0, . . . , Lr{Ur, x) = 0} 7^ 

for uq, . . . ,Ur G C"^^. This is a multihoniogeneous polynomial of degree D in each group of variables 
Ui. Thus we can directly read the dimension and the degree of V from JFy In case V is closed in P", 
its Chow form completely characterizes it, and it is possible to derive a complete set of equations 
for V from Ty. 



The main result of this work is that the computation of the Chow forms of all the equidimensional 
components of a quasi-projective variety defined by means of a given set of polynomials, has a 
polynomial complexity in terms of the number of variables, the degree, and also the length and the 
geometric degree of the input system. 

The complexity of all known general algorithms in algebraic geometry is (at least) exponential in the 
worst-case when the considered input parameters are just the number of variables and the degree 
of the input system, and there is strong evidence that this exponential behavior is unavoidable (see 
[^2| ). However, it has been observed that there are many particular instances which are much more 
tractable than the general case. This has motivated the introduction of parameters associated to the 
input system that identify these particular cases, and the design of algorithms whose performance 
is correlated to these parameters. 

In this spirit, the straight-line program representation of polynomials (sip for short) was introduced 
in the polynomial equation solving frame as an alternative data structure (see e.g. |2^, |2^, [25|) 
and its length is now considered as a meaningful parameter measuring the input (see Section |l.2| 
below for the definition of these notions and Q, 1 30 , for a broader background). 



Soon afterwards, the notion of geometric degree of the input polynomial system appeared naturally 
as another parameter classifying tractable problems. Roughly speaking, this is a parameter which 
measures the degree of the varieties successively cut out by the input polynomials (we refer to Section 



3.4 below for its precise definition). For n polynomials of degree d, the geometric degree is always 



bounded by the Bezout number d". However, there are many situations in which the geometric 



degree is much smaller than this Bezout number (see [39, Section 4] for a particular example or |3^, 
Prop. 2.12] for an analysis of the sparse case). 

In and |Q, J. Hcintz, M. Giusti, L.M. Pardo and their collaborators succeeded in classifying 
the tractability of polynomial equation solving in the zero- dimensional case in terms of the length 
and the geometric degree of the input system. They presented an algorithm whose complexity is 
polynomial in the number of variables, the degree, the length and the geometric degree of the input 
system. Our main theorem can be seen as an extension of their result to arbitrary dimensions: 

Theorem 1 Let fi, . . . , fs, g £ Q[xo, ■ ■ ■ , Xn] be homogeneous polynomials of degree bounded by d 
encoded by straight-line programs of length bounded by L. Set V := V{fi, . . . , fs) \ V{g) C P" for the 
quasi-projective variety {/i — 0, . . . , fs = 0, g ^ 0} and let V = Vq U ■ ■ ■ UVn be its minimal equidi- 
mensional decomposition. Set S := 6{fi, . . . , fs', g) for the geometric degree of the input polynomial 
system. 

Then there is a bounded probability algorithm which computes (sip's for) the Chow forms J-'y^ , . . . ,J-y^ 
within (expected) complexity s{ndS)^^^'^ L. Its worst case complexity is s{nd"')^^^^ L. 
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Let us precise the formal frame for our computations. The support of our algorithms is the model 
of bounded error probability Blum-Shuh-Smale machine (BSS machine for short) over Q: our algo- 
rithms are probabilistic BSS machines that manipulate sip's. A probabilistic BSS machine is the 
algebraic analogue of the classical notion of probabilistic Turing machine, in which the bit operations 
are replaced by the arithmetic operations {+, — , •, /} of Q. It allows to implement uniform proce- 
dures while "programming" using the basic operations of Q. This model is well suited to control the 
algebraic complexity — that is the number of arithmetic operations — performed by the algorithm. 
By hounded error probability we mean that the error probability of the machine is uniformly bounded 
by 1/4. For us, the (expected) complexity is then the expectation of the complexity seen as a random 
variable, and not its worst-case. The choice of the constant 1/4 as error probability is not restrictive: 
for any given g N we can easily modify our machine (by running it 0{\ogN) times) so that the 



final error probability is bounded by 1/N (see Proposition 1.5 and Corollary |l.6| below) 



We refer to Section 1.2 for a detailed description and discussion of the data structure and computa- 
tional model. 

We note that in our situation, it is unavoidable to consider a data structure different from the dense 
representation if we look for a polynomial time algorithm. If V is an equidimensional variety of 
dimension r defined as the zero set of a family of polynomials of degrees bounded by d, its Chow 
form !Fv is a polynomial of degree {r+1) deg V in {r+1) (n+1) variables. So the dense representation 
of J^v (i-e. the vector of its coefficients) has 

(r + l)(n + l) + (r + l) degy^^ ^ (deg F)(''+i)("+i) 



{r + l){n + l) J - ((r-|-l)(n-|-l))! 

entries, and hence it is not polynomial in degV (which in the worst case equals d"'~^). In fact. 



Corollary 2.12 below shows that in the above situation the sip representation of J-y has length 



L{Tv) = (nddegl/)<^(i)L 

The previous algorithms for the computation of the Chow forms of the equidimensional components 
of a positive-dimensional variety (|3^, [||, @) have an essentially worse complexity than ours, 
with the exception of the one due to G. Jeronimo, S. Puddu and J. Sabia ([H)), which computes 
the Chow form of the component of maximal dimension of an algebraic variety within complexity 
(s d^)^^^^ . Here, not only we compute the Chow form of all of the equidimensional components but 
we also replace the Bezout number d" by d6, where 6 denotes the geometric degree. 
Note that our algorithm also provides an effective geometric equidimensional decomposition, as each 
equidimensional component is characterized by its Chow form. 

We can easily derive equations or a geometric resolution of a variety from its Chow form (see |3.3| 
below). Hence our result contains and improves — maybe up to a polynomial — all the previous 
symbohc general algorithms for this task |ll,|2|,0, [|4|,ill,il- 

Of these, the only ones whose complexity is comparable to ours are the algorithms of Jeronimo-Sabia 
Q — which computes equations for each equidimensional component — and G. Lecerf Q, Q — 
which computes geometric resolutions. With respect to the first one, we improve the complexity by 
replacing the quantity d" by the geometric degree 6. With respect to the second one, we substantially 
improve both the probability and the uniformity aspects. Thus our work might also be seen as a 
unification and clarification of these equidimensional decomposition algorithms. 

On the contrary, it is by no means evident how to compute the Chow form from a geometric resolution 
within an admissible complexity. The difficulty lies in that the involved morphism is not finite but 



just dominant (see Remark 2.11 ). We exhibit a deter ministic algorithm which performs this task 



within complexity (sn degV)^ ' L (Main Lemma 2.3). In fact this is the key in our algorithm, and 



probably the main technical contribution of the present work. 

As a first application of our algorithm, we compute a subclass of sparse resultants. 

The classical resultant ReSn^d of a system of rt-|-l generic homogeneous polynomials in n+1 variables 
is a polynomial in the indeterminate coefficients of the polynomials, that characterizes for which 
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coefficients the system has a non-trivial solution. We show that RcSn.d can be deterministically 
computed within complexity (nd")'-'^^) (Corollary 4.1) 



The sparse resultant Res^ — the basic object of sparse elimination theory — has extensively been 
used as a tool for the resolution of polynomial equation systems (see for instance JTsI). 
Several effective procedures were proposed to compute it (see e.g. |5^, Recently, C. 

D'Andrea has obtained an explicit determinantal formula which extends Macaulay's formula to the 
sparse case (JTsI). 

From the algorithmic point of view, the basic point of sparse elimination theory is that computations 
should be substantially faster when the input polynomials are sparse (in the sense that their Newton 
polytopes are restricted). Basically, the parameters which control the sparsity are the number of 
variables n and the normalized volume Vol(^) of the convex hull of the set A of exponents. None 
of the previous algorithms computing sparse resultants is completely satisfactory, as their predicted 
complexity is exponential in all or some of these parameters (see Cor. 12.8]). 

The precise definition of Res^ is as follows: 

Let A = {ao, . . . , aAr} C be a finite set of integer vectors. We assume here that Z" is generated 
by the differences of elements in A. For i = 0, . . . , n let C/j be a group of variables indexed by the 
elements of A, and set 

aeA 

for the generic Laurent polynomial with support equal to A. Let W_a C (P^)"+i x (C*)" be the 
incidence variety of Fq, . . . , F„ in (C*)", that is 

and let tt : (P^)"+i x (C*)" ^ (pJV)"+i be the canonical projection. 

The sparse A-resultant Res^ is defined as the unique — up to a sign — irreducible polynomial in 
Z[[/o, ■ ■ ■ : Un] which defines n{WA), which is an irreducible variety of codimension 1 (see Chapter 
8, Prop.-Defn. 1.1] ). It coincides with the Chow form of the toric variety associated to the input set 



A (see Section 1.1 below). 

We show that the computation of Res^ in case A C (No)" and A contains 0, ei, . . . , e„ — the ve rtices 
of the standard simplex of R" — is an instance of our main algorithm (see Subsection 4.1.2| ). We 
thus obtain: 

Corollary 2 Let A C (No)" be a finite set which contains {0, ei, . . . , e,i}. Then there is a bounded 
probability algorithm which computes (a sip for) the A-resultant Res^ within (expected) complexity 
(n Vol(,4))'^^^^ . Its worst case complexity (nd")'^^^-', where d := max{|a| ; a G A}. 

Hence our result represents a significant improvement in the theoretical complexity of computing 
the ^-resultant as we show it is polynomial in n and Vol(^). We remark that to achieve this result, 
we had to abandon all matrix formulations. In fact, this polynomial behavior of the complexity 
is out of reach of the known matrix formulations, as in all of them the involved matrices have an 
exponential size. 

It would be desirable to extend this algorithm in order to compute general mixed resultant, a point 
will be a subject of our future research. 

As a further application, we compute the unique solution of a generic over-determined system over 
an equidimcnsional variety 1^ C P": 

Set r := dim V^. Let u = [uq, . . . , Ur) G (C"+i)''+i and set £i :— Li{ui, x) — Uio xq + ■ ■ ■ + Uin Xn for 
< i < r. 

There is a non-empty open subset of the set of coefficients u satisfying the condition that the linear 
forms £o, . . . ,ir have at least a common root in V, for which the common root is unique. It turns 
out that this unique solution ^(u) is a rational function of the coefficient u, and this fimction can be 
easily computed using the Chow form J-'y- We can successfully apply our algorithm to that situation 



see Section 4.2 for the details). 



4 



The results presented above suggest the following new research line: 

It is usual to associate to a family of s polynomials of degrees di , . . . , ds in n variables such that 
di > ■ ■ ■ > ds, the Bezout number D :— di ■ ■ ■ dmin{n,s}- The main property of this Bezout number 
is that it bounds the geometric degree of the variety defined by the input polynomials. 
We claim that the precise definition of such a Bezout number D should depend intimately on the 
encoding of the input polynomials: for polynomials of degree d encoded in dense representation, 
D := d" is a good choice, while for sparse polynomials with support in A, D := Vol(^) seems to be 
the right notion of Bezout number, as this quantity also controls the degree of the variety. 
This digression is motivated by the following crucial observation: in our computation of the resultant, 
the length of the input L together with the associated Bezout number D and the number of variables 
n controls the complexity. In case of dense representation of the input and d > 2, we have L = 
0{nd") and D = d" while for the sparse representation, we have L > 1 and D = Vol(^). In 
both cases, the complexity of computing the resultant is {nD)^'^^^ L. We think that the optimal 
complexity estimate should in fact be linear in D as well. On the other hand, it is not clear what 
should be the exact dependence on n. In the linear case, that is for n + 1 dense linear forms, we 
have L — Oin?) and D = 1, and the resultant equals the determinant, which is conjectured — but 
still not proved — to be computable in 0{n^). 

For the more general problem of computing the Chow form of polynomials encoded by sip's, the 
notion of geometric degree 5 plays the role of the Bezout number D, and our main result states 
that the complexity of computing the Chow form is s (ndD)'^^^^ L. Again, we think that the optimal 
estimate should be linear in Z?, and that the number of polynomials s and the maximum degree d 
should not occur. 

In an even more general framework, we conjecture that the computation of (a sip representation 
of) any geometric object associated to a family of polynomials in n variables represented in a given 
encoding, with associated Bezout number D and associated length of the input L, should be linear 
in both D and L, and (possibly) quadratic in n. 

Now we briefly sketch our main algorithm. 

In a first step, we prepare the input data: We take n + 1 random linear combinations of the input 
polynomials so that — with high probability — these new polynomials define the same variety V and 
behave properly with respect to the dimensions and radicality of certain ideals they successively 
define. We also take a random change of variables to ensure good conditions for the considered 
projections. 

After this preparatory step, we compute recursively the Chow forms of the components of a non- 
minimal equidimensional decomposition of V. For < r < n — 1, the algorithm deals with an 
equidimensional subvariety of the variety defined by the first n — r polynomials. 
The recursive step is as follows: 

From a geometric resolution of a zero-dimensional fiber of W^_^_i, we compute the Chow form of 
the variety obtained by intersecting W^_|_i with the set of zeroes of the next polynomial. From this 
Chow form, we obtain the Chow form of an equidimensional variety of dimension r which is a subset 
of V and contains the equidimensional component of dimension r of V together with a geometric 
resolution of the zero-dimensional fiber of that is needed for the next recursive step. 
The recursion yields the Chow forms of the components of a non-minimal equidimensional decom- 
position of V. The required minimality of the equidimensional decomposition imposes a third step 
in which we remove the spurious components. 

Finally, a word with respect to practical implementations: there is a Magma package called Kronecker 
written by G. Lecerf (see ^^) which implements — with remarkable success — the algorithm for 
polynomial equation solving in |Q, |Q (in the version of ^^). Our algorithm is closely related to 
this one and so it seems possible to implement it using this package as support. However, the lack of 
a suitable tool for manipulating sip 's is still a considerable obstacle for a successful implementation 
which would refiect the good theoretical behavior of our algorithm. 

The outline of the paper is the following: 
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In Section 1 we recall the definition and basic properties of the Chow form, and we precise the data 
structure and the computational model. We also describe some basic subroutines that we need in 
the sequel, and we estimate their complexities. 

In Section 2 we present a deterministic algorithm for the computation of the Chow form of an 
equidimensional variety from a particular zero-dimcnsional fiber, provided some genericity conditions 
are fulfilled. 

In Section 3 we describe the algorithm underlying Theorem 1, and we estimate its complexity. First 
we establish the relationship between geometric resolutions and Chow forms, and then we present 
subroutines for computing Chow forms of intersections and of components outside hypersurfaces. 
Combined with the algorithm in Section 2, this yields the desired algorithm. 

In Section 4, we apply the main algorithm to the computation of a sparse resultants, and to the 
resolution of generic over-determined equation systems. 

Acknowledgments. The authors wish to thank Agnes Szanto and Gregoire Lecerf for several 
helpful discussions, and to Joos Heintz for sharing with us his ideas concerning Bezout numbers 
and complexity. We also thank Mike Shub for suggesting us the problem of solving generic over- 
determined systems. 

1 Preliminaries 

Throughout this paper Q denotes the field of rational numbers, Z the ring of rational integers, R 
the field of real numbers, and C the field of complex numbers. We denote by N the set of positive 
rational integers, and we also denote by Nq the set of non-negative integers. 

We denote by A" and P" the n-dimensional affine space and projective space over C, respectively, 
equipped with the Zarisky topology definable over C. A quasi-projective variety V is an open set of 
a closed projective (non necessarily irreducible) variety. We denote by V C P" the projective closure 
of V, that is the minimal closed projective variety which contains it. 

If /i, . . . , /s are polynomials in Qfxo, . . . , a;„], V{fi, . . . , fs) will denote the set of common zeros of 
fi, . . . , fs in P". This notation will also be used in the affine case. 

Let 1^ be a quasi-projective variety and let g e Q[a;o, . . . , x„] be a homogeneous polynomial. Then 
we denote by Vg the basic open set V \ V{g) of V. 

We adopt the usual notion of degree of an irreducible projective variety. The degree of an arbitrary 
projective variety is here the sum of the degrees of its irreducible components. If the variety is 
quasi-projective, its degree is defined as the degree of its projective closure. 

We only consider polynomials and rational functions with coefficients in Q and varieties defined by 
polynomials with coefficients in Q unless otherwise explicitly stated. 

The determinant of a square matrix M is denoted by |M|. 

1.1 The Chow form of a quasi-projective variety 

We gather in this subsection some definitions and basic facts about Chow forms. For a more detailed 
account we refer to |4^, Section 1.6.5], Chapter 3] and Q. 

First we define the notion of Chow form of an equidimensional quasi-projective variety: 
Let C P" be an equidimensional quasi-projective variety of dimension r. 

For i = 0, ... ,r let Ui = {Uio, Un, . . . , Um) be a group of n -I- 1 variables and set U := {Uq, . . . , [/„). 
Then set 

Li := [/jo xo H h U^n Xn eQ[U] [x] 
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for the associated generic linear form, where x denotes the group of variables [xq, . . . Let 

<i>v = {{uo, ■■■,ur;0& X P" ; e e ^, ^0(^0, = 0, ■ ■ ■ , Lr{Ur,0 = 0} C (P")''+i X P" 

be the incidence variety of these linear forms in V, and let tt : (P")''+i x P" ^ ^pny+i -^^ ^-j^^ 
projection (u,^) u. 



Lemma 1.1 Under the previous assumptions and notations, 7r($y) = 7r(<i>y). 

Proof.- Let V = Uc C be the irreducible decomposition of V . From the definition above we deduce 
that $y = Uc $c and so 7r($y) = Uc 7r($c)- 

We also have that V = Uc is the irreducible decomposition of V . Then tt{^y) — Uc 7r($^) and 
so, without loss of generality, we can assume that V is irreducible. 

The map V defined by (u,^) (, makes a fiber bundle over V with fiber (p"-i)'+i. 

Then <i>y is an irreducible variety of codimension n + 1, and the same is true for $-^7. 
As is a closed set, $y C ^y. These are irreducible projective varieties of the same dimension 
and, therefore, they are equal. The fact that tt is a closed map implies that 7t{^v) — T^i^Y)- '-' 

Then 7r($y) C (P"-)'"+i is a closed hypersurface p. 66]. We define a Chow form of V as any 
squarefree defining equation J^v & Q[Uo, ■ ■ ■ , Ur] of the Zariski closure 7r(<i>y) C (P")''+^. Note that 
the Chow form of an equidimensional variety is uniquely determined up to a scalar factor. We extend 
this to dimension —1 defining a Chow form of the empty variety as any non-zero constant in Q. 
This definition extends the usual notion of Chow form of closed projective equidimensional varieties. 
In fact. Lemma states that a Chow form of an equidimensional quasi-projective variety is a Chow 
form of its projective closure. 

From this definition, we see that any Chow form of V characterizes the sets of over-determined linear 
systems over the variety V which intersect it: for uq, . . . , Ur £ C"+^ we have 

JPy(wo,.-.,Ur) = ^ Fn {io(uo,a;) = 0} n • •• n = 0} 7^ 0. 

A Chow form J-y is a multihomogeneous polynomial of degree deg V in each group of variables 
Ui (0 < i < r). The variety V is uniquely determined by a Chow form of V i W^, P . 66]). Moreover, 
it is possible to derive equations for the variety V from a Chow form of V ([[19|, Chapter 3, Cor. 
2.6]). 

In case V is irreducible, J-y is an irreducible polynomial and, in the general case, a Chow form of 
V is the product of Chow forms of its irreducible components. 

Following we avoid the indeterminacy of J-'v by fixing one of its coefficients under the following 
assumption on the equidimensional quasi-projective variety V: 

Assumption 1.2 If dim V ~ 0, we assume V C {xq ^ 0}. //dimy = ?■ > 0, we assume that the 
projection iry : V P'' defined by x ^ (xq : • • • : Xr) verifies 4I^TTy^{{l : : • • • : 0)) = deg V . 

This assumption implies that V n {xi = 0} H • • • H {xr = 0} is a 0-dimensional variety lying in the 
affine chart {xq 7^ 0}. In particular V has no components contained in the hyperplane {xq — 0}. 
We also note that, in case is a closed affine variety, the hypothesis #(V^ H {xi = 0} n • • • n {xj. — 
0}) — deg implies that the map Try : V ^ hJ" defined by x 1— > {xi, . . . , Xr) is finite; that is, the 
variables xi, . . . ,Xr are in Noether normal position with respect to V Lemma 2.14]). 
Set Ci for the {i -f l)-vector of the canonical basis of Q"+^ and D := deg V. Then, under Assumption 



1.2, !Fv{so, . . . , Gr) — that is, the coefficient of the monomial U^q ■ ■ ■ U^j. — is non-zero. Then, we 



define the (normalized) Chow form Chy of V by fixing the choice of jFy through the condition 

Chvieo, . . . , gj.) = 1. 



Note that if V satisfies Assumption |L2| then each of its irreducible components also does. Therefore, 
the normalized Chow form of V equals the product of the normalized Chow forms of its irreducible 
components. The normalized Chow form of the empty variety equals the polynomial 1. 
Here there are some examples of Chow forms: 
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• In case dim V ^ we have 

Furthermore, if V satisfies Assumption |L.2| , Chy is equal to the above expression provided we 
choose homogeneous coordinates of the type ^ := (1 : f') G P" for each point in V. 



In case ^ is a hypersurface of degree d, then V = V{F) C P" where F e Q[xo, . . . ,Xn] 
is a squarefree homogeneous polynomial of degree d. We consider the n x (n + l)-matrix 
M := {Uij)i<i<n, and, for < j = Q < n, we set Mj for the maximal minor obtained by 

0<j<n 

deleting its (j + 1) column. Then 

Tv = F(Afo, -Ml, . . . , (-1)" Mn) e Q[Uo, ...,[/„]• 



In this case, Assumption 1.2 is equivalent to the fact that / :— 0, . . . , 0, is a squarefree 
polynomial of degree d in t. Therefore, Chy is equal to the above expression if we choose F 
such that the coefficient of the monomial x'^ is 1. 

• The sparse resultant provides an important family of examples: let A = {ao, . . . , ajv} C Z" be 
a finite set of integer vectors, such that the differences of elements in A generate Z". Consider 
the map 

V?^ : (C*)" P^ , ■■■■■■■ C")- 

This is always well-defined as ^ (1 < i < n) for all ^ G (C*)". The Zariski closure of the 
image of this map X_a := (^^((C*)") C P^ is the toric variety associated to the set A. This is 
an irreducible variety of dimension n and degree Vol(,4) (the normalized volume of the convex 
huh of A). 

The ^-resultant equals the Chow form of this variety Chapter 8, Prop. 2.1], that is: 

I'xa = R-es^. 

We refer to and to Chapter 7] for a broader background on sparse resultants and toric 
varieties. 



1.2 Data and algorithms structure 

First we specify our data structure: 

The objects we deal with are polynomials with rational coefficients. The data structure we adopt to 
represent them concretely is the straight-line program encoding (sip's for short). The input, output 
and intermediate objects computed by our algorithms are polynomials codified through sip's. We 
insist on the fact that in the present work the crucial feature of sip's is their role as data structures, 
rather than the more traditional functional role as programs without branchings for the evaluation 
of polynomials at given points. 

For the standard terminology of sip's, see Defn. 4.2]. In this paper all sip's are defined over Q, 
without divisions and expecting the variables xi, . . . ,x„ as input. 

For purpose of completeness we restate the definition in our particular case: 

Let n G N. We denote by {+, — , •} the addition, substraction and multiplication in the Q-algebra 
Q[xi, . . . ,Xn]- We consider apart the addition and multiplication by scalars, that is for A G Q and 
/ G Q[xi, . . . , a;„] we set A^(/) := / + A and A"^(/) := A • /. We denote by Q'^ and Q"^ the set of all 
scalar additions and multiplications for A G Q, respectively. 

We set r2„ := U U {+, — , •} and denote by ar(u;) the arity of an operation u G fi„: that is 1 
if it is a scalar operation and 2 if it is a binary one. 
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A straight-line program 7 (over Q and expecting xi, . . . ,Xn as input) is a sequence 7 :— (71, . . . , 7^) 
of instructions 

_ f (wi; fcii) if ar(a;i) = 1 
~ \ {uj,; kii,ka) if ar(wi) = 2 

where each G fin is an operation and for every j, /cij £ Z satisfies — n + 1 < fc^- < i — 1 and 
represents a choice of a previous index. The number of instructions L in 7 is called the length of 7 
and is denoted by L{j). This is, in the standard terminology, the complexity induced by the cost 
function which charges 1 to each operation in Jl„ (see |7[ Defn. 4.7]). 

Given a sip 7 ~ (71, . . . , 7^), its result sequence (/-„+i, . . . , /l) is classically defined as 

Here we make a slight modification with respect to the standard terminology. According to the data 
structure role played by sip's we consider only the final result of the sip 7, that is the final polynomial 
Jl G Q[xi, . . . , Xn]- We call it the result of 7. Here is the precise definition: Let Tq[xi, . . . , Xn] denote 
the set of sip's over Q expecting xi, . . . , a;„ as input. Then there is a well-defined surjective function 

Eval : rQ[a::i, . . . ,a;„] (Q[a;i, . . . ,x„] , 1 ^ fb where L := L{-f). 

In this way each sip defines precisely one polynomial (and not a finite set). We say that 7 G 
rQ[xi, . . . , Xn] encodes f € Q[xi, . . . , Xn] if / is the result of 7. Given a polynomial / G Q[a;i, . . . , x„] 
we define its length L{f) as the minimal length of a sip which encodes /. (We always have deg/ < 
2^(/).) For a finite set V C Q[xi, . . . , Xn] we define naively its length as L{V) :— '^f^-p L{f). 

From the dense representation Oq, of a polynomial / G Q[xi, . . . , Xn] we obtain straight- 
forward a sip for /. First, it is easy to show inductively that for any r G N, there is a sip of length 
bounded by ("^'^) whose result sequence is the set of all monomials x" of degree |a| < r. This is 
due to the fact that once one has a list of all monomials of degree bounded by r — 1, each one of the 
homogeneous monomials of degree r is simply obtained from one of the list multiplying by 
a single variable. Now set d :— deg/. We multiply all monomials of degree bounded by d by their 
coefficients and add them up, that is we add 2 ("J'') instructions to the sip, in order to obtain a sip 
which encodes /. Hence 

i^(/)<3(^" + ^) <3(d+l)". 

We call this particular sip the dense sip of /. The previous computation shows that in the worst 
case, the length L{f) of a polynomial / of degree d is essentially its number of monomials. 

We can operate with the data structure sip, extending directly the operations in f2„: for instance 
for * £ {+, — , •}, given two sip's 7, (5 6 rQ[xi, . . . , Xn] we obtain the new sip 

7 *(5 := *(7,(5) := (71, . . . , 7i(^), ^i, . . . ,(5^(5), (*;i(7),i(7) + L{S))), 

where the choice of previous indexes for S are suitable modified. This sip obviously encodes the * of 
the two polynomials encoded by 7 and 6, and its length is £(7) + L{S) + I. 

More generally, for 7 G rQ[2/i, . . . , and Si, . . . ,Srn G rQ[a;i, . . . , x„], we can define the composition 
sip 7 o (5 := 7 o ((5i, . . . ,Sm) G rQ[a;i, . . . ,x„]. We have 

L(7 o d) = L{6i) + ■■■ + LiS,n) + i(7)- 

This operation is compatible with the map Eval, that is Eval(7 o 5) — Eval(7) o Eval(5). Hence 

for / G Q[yi, . . . ,y„i] and 51, . . . ,5„ G Q[xi, . . . ,a;„] we have that £(/(gi, . . . ,gn)) < L{gi) H h 

Lign) + L{f). 

Now we specify the computational model that produces and manipulates our data structure: it is 
the Blum- S hub -S male machine over Q (BSS machine for short), which captures the informal notion 
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of uniform algorithm over K. We refer to |Q, Chapters 3 and 4] for the definition, terminology and 
basic properties. However, there are again some slight modifications in our definition (restrictions 
on the operations -only over rational numbers- and the branchs -only equality of numbers to zero), 
and we restate it for purpose of completeness: 

We recall that a BSS machine M over Q has five types of nodes: input, computation, branch, shift, 

and output. 

Set 

:= y 

n>0 

for the disjoint union of the n-dimensional spaces Q", i.e. the natural space to represent problem 
instances of arbitrarily large dimension. For a G C Q°° we call £ the size of a, and we denote it 
by size (a). 

On the other hand, let 

Qoo := Q 

be the bi-infinite direct sum space over Q. The elements h G Qoo arc of the form 

b:= {... ,6_2,&-i,feo ■ 6i,&2, • ■ • ): 

where bi — Q for \i\ 3> 0. The dot between 6o and bi is a distinguished marker which allows to 
visuaHze the position of the coordinates of b. 

Now we define the computation maps: For each uj G {+, — , ■, /} and i, j. A: G N there is a map 
Qoo ^ Qoo , b^ {. . . ,bk-2,bk-i,i^[bi,bj),bk+i,. ■ ■ ) 

(observe that unlike in the case of our data structure, here we also allow division). Also for each 
A° G Q^ or A'" G Q™ and i. A; G N there is in an analogous way a map 

Qoo ^ Qoo , b^[...,bk-2M~U>^''^'"\h),bu+l,...). 

These are our computation nodes. 

The only branch node we allow is the one which checks the equality 5i = 0. In other words its 
associated map is Qoo ^ Q , 6 i— > 6i. 

The shift nodes are of two types: shifting to the left or to the right, associated with the two maps 
Qoo ^ Qoo , b^ (7i{b)i^bi+i or (7r{b)i^bi^i. 

The machine M. over Q is a finite connected directed graph containing these five types of nodes 
(input, computation, branch, shift and output). The space Q°° is both the input space Im and 
the output space Omi and Qoo is the state space Sm, that is the "working" space of the machine. 
The dimension Km of the machine is the maximum dimension of the computation maps, which, 
under our assumptions, coincides with the maximum of the natural numbers i,j and k involved in 
the computations. 

We are interested in the algebraic complexity of these machines. We assume that the cost of each 
computation and branch node is 1, and that of the shift nodes is 0. In other words, we assume 
that M can do arithmetic operations and equality questions in Q at unit cost, and shifts at no cost. 
Hence the complexity CMifl) of the machine M. on an input a is just the number of computation 
and branch nodes of the graph, from input to output. Note that each computation node performs 
exactly one operation, and hence C7K(a) is just the number of basic arithmetic operations and 
equality questions performed on the input a. 

Observe that any sip 7 G rQ[a::i, . . . , is an example of a (finite-dimensional) BSS machine 
without branchs or shift nodes for computing / :— Eval(7) G Q[a:i, . . . , x„] at any input point 
a G Q". The dimension of this machine is n + ^(7) and its complexity is ^(7). 
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Given £ 6 N we consider the complexity Cm (^) of the machine on inputs of size bounded by £, that 
is 

Cm{^) ■= sup{Cx(a) ; size(a) < £}. 



This computational model is a natural algebraic analogue of the notion of Turing machine. In 
particular it provides us a support for the implementation of uniform procedures. 
The difference with the Turing model is that we replace bit operations by arithmetic ones, and that 
we do not count the cost of the shifts operators. The reason not to count the shifts is not fundamental, 
and is just that the available literature on the algebraic complexity of basic polynomial manipulation 
only takes into account the number of arithmetic operations and comparisons (see Subsection 1.2, 
see also 0). 

Since all the involved computations arc done over the rational field, the machine Ai can be effectively 
transformed into a classical Turing machine. However our complexity counting docs not provide any 
reasonable control on the complexity of the resulting Turing machine. 



As we have already anticipated, our algorithms are BSS machines that manipulate sip's. A machine 
A4 receives as input a finite family of sip's 7 G r(j[a;i, . . . and gives back a finite family of 

slp's7W(7) erQ[2/i,...,j/,„]*-^'. 

A finite family of sip's ^ GTqIxi, . . . , x„]*^ can be easily codified as an input element in Im = Q°°, 
in fact 7 can be identified with a vector in the following way: 

The first coordinate is for the dimension n, that is the number of variables. Then each instruction 
of the first sip 71 is codified as a triple: we enumerate the operations in fin with numbers from 
2 to 6, 2 and 3 corresponding to the operations in Q° and Q™, and 4 to 6 to +, — and •. For 
operations in <Q^ U we write the operation number in the first coordinate, the corresponding 
coefficient in the second one, and the position to which it applies in the third one. The binary 
operations are codified in a similar way, by writing first the operation number, and then the position 
of the two elements to which it applies. The positions are numerated from 1 — n to i(7) according 
to the definition of result sequence. For instance the vector (2, (3, 5, —1), (4, 0, 1), (6, 2, 2)) codifies 
the sip xi, a;2, 5x1, X2 + 5xi, {x2 + Sxi)^. The instruction to separate two consecutive sip's is an 
empty cell, that is a 0. The second sip 72 is now codified exactly as the first one. Therefore, 
7 := (71, . . . ,7Af) e rQ[xi, . . . ,Xn]^ is codified as a vector in Q°°, in fact in Q^^+^^Ct) since we 
need to add M — 1 "0" to separate two consecutive sip's. 

The machine 7W manipulates this input, the finite family of sip's 7 G rQ[xi, . . . ,x„]^^: it operates 
these sip's and gives as the output an element of Om corresponding to a finite family of sip's in 
rQ[yi, . . . ,?/m]*^ ■ As we've just seen, the input and output size is (essentially) the length of each 
of these families. Thus, we speak here of a finite family of sip's 7 as the input of M. and we simply 
denote by Aiij) its output in rQ[2/i, . . . , j/m]*^ • 

Remark 1.3 Let 7 G rQ[a;i, . . . , Xn]^'^ be the input sip family of a BSS machine M and let A^(7) G 
rQ[yi, . . . , be its output. Then 

L{M{^))<iL{^) + CM{l)- 

Proof.- As we do not know how the machine M operates on 7, the only bound for i(A^(7)) is the 
number of operations labelled from 2 to 6 of the representation of M{'^) in A^, which is bounded 
by the number of non-zero cells of this representation minus 1 (since the first cell of the output 
corresponds to the number of variables m of the output) . This is bounded by 1 + 3 £(7) + Cm (7) — 1 j 
that is, the size of 7 as an input of (excepting the M — 1 zero-cells separating different input 
sip's) plus the number of computation nodes Cm{i) minus 1. □ 



Our main algorithms are in fact probabilistic. For this reason we implement them in the model of 
probabilistic BSS machine over Q Section 17.1]. This is a BSS machine Ai with an additional 
kind of node, called probabilistic. These are nodes that have two next nodes and no associated map 
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and that "flip coins" , that is when a computation reaches a probabihstic node, it randomly chooses 
the next node between the two possible ones with probability 1/2 for each of them. 

In this probabilistic setting, each run — on the same given input 7 — of the machine Ai may lead 
to a different path computation. In our case, for any given input, the number of probabilistic nodes 
traversed is finite, that is, the number of possible paths is finite. We treat the probalistic nodes as 
branches and charge cost 1 for each of them. 

As every path V oi M. corresponds to a BSS machine of complexity C-p{'y), the algebraic complexity 
Cm{i) of the machine Ai on the input 7 turns out to be a random variable, with finite sample set. 
Moreover, again in our context, every path is finite: it may happen that a branch ends in an error 
message but in any case the complexity of any path is bounded. Thus the random variable Cm (7) 
satisfies 

Prob {Cm{i) = C):=Y, Prob(P; V path such that Cp(7) = C). 

We are interested in the worst case complexity C^^{'y), the maximum complexity of the paths of 
on 7, and the expected complexity E^ij), defined as the (finite) expectation of this random 
variable, that is 

Em{i) := E{CMil)) = E • Prob(C>,(7) = C). 

Observe that C'^^i^) > Em{i) always holds. 

As before, we also consider the function Em : N — *■ N: 

EM{i) ■■= sup {Em{i) ; n, M e N, 7 e rQ[a;i, . . . , a;„]^ and M + 3 1.(7) < l). 

We define now the error probability of the machine on a given input. Again, there is here a modifi- 
cation with respect to traditional probabilistic BSS machines. Keeping in mind that for any run of 

the probabilistic machine M. on the input 7 G r(g)[.Ti, . . . , .t„]'*^, its output (independently from the 
path randomly taken) encodes a finite family of polynomials / := (/i, . . . , /m) G Q[yi, . . . , ym]^ 
we define: 

Definition 1.4 (Bounded Probability Algorithm) Given 7 G rQ,[xi, . . . ,Xn]^ and given a set of 
polynomials f := (/i, . . . , /m') G Q[yi, ■ ■ ■ , ym]^ , the error probability ex(7, /) that M computes 
f on the given input 7 is the probability that the output of A4 onj does not encode f ; that is the prob- 
ability that the computation finishes with an error message, or that it outputs d G Tqlyi,. . . ,ym\^ 
which does not encode f. 

We say that M. computes / if ex (7, /) < 1/4. As this happens at most for one f, when it happens, 

we set CM (7) := il, f) ■ 

When the latter happens for every input, we say that M is a bounded probability machine for 
polynomial sip's, and we speak 0/ bounded probability algorithm. 

Observe that our probabilistic machine is a little unusual since in fact, as different sip's may encode 
the same polynomial, the polynomial / computed by the probabilistic machine M corresponds to 
an equivalence class of outputs rather than a single one. In this paper, all machines are bounded 
probability machines for polynomial sip's in the sense of this definition. 

In our setting, probability is introduced by choosing a random element with equidistributed prob- 
ability in a set [0,^)" := {0,...,£— l}" for given natural numbers £ and n. Since probabilistic 
machines flip coins to decide binary digits, each of these random choices can be easily simulated 
with a machine with complexity 0{n \log£]), where here and in the sequel, log denotes logarithm 
in base 2. This machine is denoted by Random(n,£). 

In this work, in each case, there is a non-zero polynomial F gQ[xi,. . . , Xn] \ {0} such that a random 
choice a is good — that is, leads to the computation of the desired output — if F{a) 0. The error 
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probability of this random choice is then estimated by the Zippcl-Schwartz's zero test ( 53 
Lemma 1]): 

dee F 

¥mh{F{a) = 0) < . 

The choice of 1/4 as a bound for the error probabihty is not restrictive and we can easily modify it in 
order to reduce the error probability as much as desired. The usual procedure is to run the machine 
many times and to declare that the polynomial family / is computed by M. if it is the output of 
more than half the times. There is a slight difficulty here, appearing from the fact that our machine 
computes sip's instead of polynomials, and two different runs may lead to different encodings of the 
same polynomials. That is why we need here to be more careful in our definition. We define it in 
the following way: 

Given the bounded probability machine M which on input 7 e rQ[a;i, . . . , a;„]*'^ computes / e 
Q[yi, . . . , Um]^ , and given s G N, the machine Ads is the machine which performs the following: 

1. A4s runs s times the machine M on the given input 7: for 1 < i < s it obtains the output sip 
family Si together with the complexity Ci of the path followed to compute Si. 

2. Then Ms chooses randomly a £ [0, M' 2'*+^ and computes S^{a), l<i<s. 

3. For j = 1 to [s/2] 

• it computes Sj{a) — Sk{a),j < k < s, and compares the results to 0. 

• if is obtained for strictly more than s/2 values of fc, it returns the polynomial family 
encoded by Sj as the output and ends. 

• if not, it goes to j + 1 

If for no j < [s/2] the algorithm obtains enough times, it outputs error and ends. 

Proposition 1.5 Given a bounded probability machine A4 which on j ^ TqIxi, . . . ,a;„]^^ computes 
f £ ^[yii • ■ • , y™]^^ <^'>T'd given s G N, the worst case complexity, the expected complexity and the 
error probability of the machine A4s on 7 verify the following bounds: 

C^T(7) = 0((™ + l)s(i(7) + C^-(7))+mlogM')+A/'Q), 

EmAi) - 0((m + l)s(L(7)+i;^(7)) + mlogAf' + M'Q), 
cmM < 2(3/4)^/2. 

Proof - Let us first describe the algebraic complexity C of a given run of the machine Cm, on 7 in 
terms of the complexities Ci of the paths followed by the machine on 7 on the i-th run. 

1. has complexity Ci + ■ ■ ■ + Cs- 

2. Producing the random choice a costs 0{m (logM' + s + 1^(7) + Ci + ■ ■ ■ + Cs)), and, from Remark 



1.3| , computing Si{a), . . . , Ssia) costs 3s L{"f) + Ci + ■ ■ ■ + Cs- 

3. As Sj{a) e Q^^ , to compute all Sj{a) — 5k{a) and compare them to costs 2 M' (2). 
Hence, the worst-case complexity of the machine Ms on 7 is 

CZrh) = 0{{m + 1) s (L(7) + CZ^-^ij)) + m logM' + M' (^^ 

while, as the complexity is an affine combination of the s independent random variables Ci, . . . , Cg 
its expectation verifies 

Em. (7) = 0{{m + 1) s {L{j) + Em{i)) + m logM' + M' (^^^ )• 
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The error is bounded by the probabihty that there is no group of more than s/2 vectors which 
coincide plus the probabihty that Sj (a) = Sk (a) but the two polynomial families encoded by Sj and 
Sk do not coincide. 

The first error is bounded by (3/4)*/^ as in Sect. 17.2, Lemma 1]. To estimate the second error 
we apply Schwartz' lemma: for 1 < « < s the output 6i codifies / G Q[j/i, . . . , j/m]*^ where the 
degree of each component is bounded by 2^^^''')+'-'*. Thus the error of one test is bounded by 
(M'23^(')')+'^')/(A/'2"+3^(t)+Ci+---+c=) < (1/2)^ As there are at most Q such independent tests, 
the total error verifies 



eA..(7)<(3/4r/'+(2) (l/2r <2(3/4r 



for s > 2. □ 

Corollary 1.6 Given a bounded probability machine M which on j ^ TqIxi, . . . , Xn]^^ computes 
f G Q[yi; • • ■ J Um]^'^ o,nd given N € N, N > A, the error probability of the machine Ais on 7 for 
s := [6 (log + 1)] is bounded by 1/N while its worst-case complexity is of order 

0((m + 1) logiV (i(7) + C^""(7)) + m log Af + M' log' N). 

Proof - As (3/4)3 < ^/2), 

eMAl) < 2(3/4)3 0°s^+i) < (3/4)3 i°g^ < (l/2)i°s^ = l/N. 

□ 



Proposition 1.6 will be used to decrease the error probability of intermediate subroutines of our 
main algorithm and keep control of the complexity in order that the error probability of the latter 
is bounded by 1/4. Observe that the length of the output sip is of the same order that the length 
of the sip obtained when running any of the repetitions of the algorithm. 

Given a bounded probability machine A4, any time we want to obtain the output of A4 for a sip 
input family 7 with error probability bounded by l/N, we run Subroutine |l| which gives a new 
probability machine A^(7; N) doing so. Any time we run M for the input family 7, we will denote 
by Complexity (A^ (7)) the complexity of doing it this time. 
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Subroutine 1: Decreasing error probability of M. 



procedure A4{'^;N) 

^ 7 is a sip input family for M. 
# e N,iV > 4. 

=ff The procedure returns the output of A4 with error probability bounded by 1/N. 

1. s := [6(logAr+ 1)1; 

2. for i from 1 to s do 

3. [S,, a) (Mij), Complexity(A^(7))); 

4. od; 

5. a := Random(m, M' 2'*+3i(T)+'^i+-+c.); 

6. (Ai,...,A,):-(<5i(a),...,(5,(a)); 

7. J ■■= 1; 

8. while j < [s/2] do 

9. k:^j + 1; 

10. i:=0; 

11. while k < s do 

12. if Aj - Afe = then 

13. t:^t+l; 

14. fi; 

15. fc:=fc + l; 

16. od; 

17. if t > s/2 then 

18. return(5j); 

19. else 

20. J := j + 1; 

21. fi; 

22. od; 

23. return( "error"); 
end. 



1.3 Complexity of basic computations 

We summarize the complexity of the basic operations on polynomials and matrices our algorithms 
rely on. We refer to for a rather complete account of the subject or to |2^ for a brief survey of 
the existing literature. 

Let R denote a commutative Q-algebra and let d e N. 

The multiplication of d x d-matrices with coefficients in R can be done with 0{d? '^^) operations 
of R and no branches (||^, Cor. 15.45]). The computation of the coefficients of the characteristic 
polynomial of a d x d-matrix — and in particular the computation of the adjoint and the determinant 
of this matrix — can be done with 0{d'^) arithmetic operations and no branches, the same bounds 
hold for the inversion of an invertible matrix ([0, H]). 
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The quantity 

M{d) := d \ogd log log d 

controls the complexity of the basic arithmetic operations (addition, multiplication, and division 
with remainder) for univariate polynomials with coefficients in R of degree bounded by d in dense 
representation: Addition of univariate polynomials can be done in d+ 1 arithmetic operations, while 
the Schonhage-Strassen polynomial multiplication algorithm takes 0{M{d)) arithmetic operations 
and has no branches ([Q, Thm. 2.13]). Division with remainder — provided the divisor is a monic 
polynomial — has also complexity 0{M{d)) and no branches (||^, Cor. 2.26]). The greatest common 
divisor can be computed through subresultants with 0{d) branches (computing the degree of the gcd 
corresponds to checking the vanishing of the determinant of submatrices of the Sylvester matrix) and 
complexity 0{d^) (solving a linear system) ((l2j, Q). Alternatively, the Knuth-Schonhage algorithm 
could be used to compute the greatest common divisor with complexity 0{AI{d) logo?) and 0{d^) 
branches (|0, Cor. 3.14]). 

Now we are going to consider some procedures involving polynomials encoded by sip's. 
First, given a sip 7 which computes / G Q[a;i, . . . , Xn] and given a G Q", we can compute /(a) G Q 
within complexity ^(7) and so, we can also check /(a) = within the same complexity. The 
derivative of the polynomial / with respect to one of its variables will be computed by means of the 
Baur-Strassen's algorithm (see |^) within complexity 0(L(7)). 

For a group of variables y := {yi, . . . ,ym) and a £ Q™, we will denote by Expand(/, y, o, c?) the 
subroutine which, given a sip 7 which encodes a multivariate polynomial /, computes as intermediate 
results sip's for the homogeneous components centered at a and of degree bounded by d of the 
polynomial / with respect to the given group of variables y: In Lemma 13], Lemma 21.25] are 
given sip's of length 0{d^ -^(7)) in which all the homogeneous components of / of degree bounded 
by d appear as intermediate computations. These procedures can be easily modified within the 
same complexity to compute the homogeneous components centered at a and up to degree d of 
a polynomial with respect to the given group of variables. In particular, if y consists of a single 
variable and a — 0, this procedure computes the coefficients of the given polynomial with respect to 

y- 

Quite frequently we use a mixed representation of /: instead of encoding it by means of a single sip, 
we consider / as a polynomial in a distinguished variable, and if d is a bound for the degree of / in 
this variable, we give a {d+ l)-uple of sip's, which encode the coefficients /o, . . . , /d of / with respect 
to the variable. The length of this mixed encoding does not essentially differ from the length of /; 
denote by L'{f) the length of the mixed encoding and by L{f) the length of /, we have: 

L{f) = 0{d + L'{f)) and L' {f) ^ 0{d^ L{f )). 

Sometimes we need to compute the exact degree of a polynomial with respect to a particular variable. 
We will call Deg(/, d) the procedure which computes the degree of the univariate polynomial / 
given by its dense representation, where d is a bound for its degree. This computation is done 
by simply comparing the coefficients of / with 0. This procedure can be adapted to obtain a 
probabilistic algorithm Deg(/i, . . . , fs, x, d; N) which computes, with error probability bounded by 
1/N, the total degrees of the polynomials fi, ■ ■ ■ , fs in the group of variables x, from sip's encoding 
/i, . . . , /s and an upper bound d for their degrees in the variables x. To do so, first we apply 
subroutine Expand(/i, x, 0, d) for 1 < i < s, to obtain the homogeneous components of fi. Then 
by choosing a random point in [0, 1, ... , sdN)" we decide probabilistically which is the component 
of greatest degree different from zero of each polynomial /i, . . . If the given polynomials are 
encoded by sip's of length boimded by L, the worst-case complexity of this procedure is of order 
0{sd^L + nlog{sdN)). 

1.4 Effective division procedures 

Here, we gather the division subroutines we will need. Basically, they compute the division of multi- 
variate polynomials and power series, and the greatest common divisor of multivariate polynomials. 
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Subroutine 2: Polynomial Division 



procedure PolynomialDivision(/, g, d, a) 

f-<9 ^ Qia^i, • ■ • , Xn] such that / divides 

# d G N an upper bound for the degree of the quotient g//, 

# a G Q" such that /(a) ^ 0. 

# The procedure returns h :~ g/ f . 
1. «:=/(a); 

3. H^g-v{a~f); 

4. (iJo, . . . , := Expand(i/, X, a, d); 

5. h:— '}2m=o ^-m] 
return(/i); 

end. 



In all cases, the objects will be multivariate polynomials encoded by sip's and power series, whose 
known graded parts will be also encoded by sip's. The proposed procedure for multivariate power 
series division is new and plays an important role in Subroutine 0, which in turn is the key step of 
our main algorithm. 

The following subroutine is the well-known Strassen's Vermeidung von Divisionen (division avoiding) 
algorithm ( |p9| ). We re-prove it briefly in order to estimate its complexity. 

Lemma 1.7 Let f,g ^ Q[xi, ■ ■ ■ ,Xn] be polynomials encoded by sip's of length bounded by L such 
that f divides g. Let d € N be such that deg{g/ f) < d, and a G Q" such that f(a) ^ 0. 
Then PolynomialDivision (Subroutine^ computes g/f within complexity 0{d^ {d + L)). 

Proof.- The quotient polynomial h := g/f G Q[ also be seen as a power series in 

Q\\x — a]]. For a :— f{a), we have 

For 



=0 



d 



i=0 

we have h = H mod [x — a)''+^. Thus, if {Hm)m<d are the homogeneous components of H centered 
at a and of degree bounded by c? , we conclude h = Ylta=o 

The stated complexity is obtained as follows: We compute the univariate polynomial v with ©(d+L) 
operations. Hence we compute H within complexity 0{d + L). We compute its homogeneous 
components in x — a up to degree d within complexity 0{d'^{L{H)}) = 0{d^{d + L)). Finally we 
obtain h as X]m=o within the same complexity bound. □ 

Observe that the same procedure can be used to compute the graded parts centered at a and of a 
certain bounded degree of the rational function 5//, even in case / does not divide g. We denote 
this subroutine by GradedParts(/, g, Z?, a), where the argument D corresponds to the bound for the 
degree of the graded parts to be computed. Its complexity is of order 0{D^{D + L)). 
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Subroutine g converts sip's with divisions computing polynomials in Q[a;i , . . . , Xn] into ordinary sip's: 
Sip's with divisions are defined as ordinary sip's, but with the set of basic operations enlarged to 
include the division, which we denote by the bar /. A further requirement is that all divisions should 
be well-defined, that is, no intermediate denominator can be zero. In general, the result of a sip with 
divisions is a rational function in (Q)(a;i, . . . , 

Observe that, given a sip with divisions 7 which encodes a rational function h, we can easily compute 
separately a numerator g and a denominator / by means of two sip's (, 77 without divisions: for 
instance, for each addition hi := hj + hk in the result sequence of 7, if hj := hji/hj2 and h^ '■— 
hki/hk2, we set gk ■= hji hk2 + hj2 hki and hk := hj2 hk2 for the corresponding result sequence in ( 
and T] respectively. We proceed analogously for the other operations in fJ^ := fin U {/}. 
We have 

Furthermore the sip's ^ and rj can be computed within complexity L(C) < 3-^(7) and L{ri) < L{'~f). 
In particular, given a £ Q", we can check if 7 is well-defined at a and, if that is the case, if h{a) — 
within complexity 0{L{j)). 

In case /i is a polynomial of degree bounded by d, the previous considerations together with Lemma 



1.7 show that we can compute a sip without divisions for h with complexity 0{d?'(d + £(7))) 



Now follows a bounded probability algorithm (in the sense of Definition 1.4) to compute the greatest 
common divisor of two multivariate polynomials encoded by sip's (for another algorithm solving this 
task see |^). Herein, GCDl(i^, G, d, e) is the subresultant algorithm which computes a greatest 
common divisor of two univariate polynomials F and G of degrees d and e respectively with coeffi- 
cients in a ring A. The output (g, Q) of GCDl is the multiple Q modulo the factor q G A — {0} of 
the monic greatest common divisor of F and G over the fraction field of A. 

Lemma 1.8 Let f,g£ Q[a;i, . . . , Xn] be polynomials of degree bounded by d encoded by sip 's of length 
bounded by L. 

Then GCD (Subroutine^ is a bounded probability algorithm which computes (a sip for) the greatest 
common divisor between f and g. Its (worst-case) complexity is of order 0{n log d -\- d'^ (d^ -\- L)). 

Proof.- For a £ Q" such that /(a) 7^ and t an additional variable, we set 

i^(^,t) :=i^/(^+a) , G{x,t):^t^gi^+a) e Q[x][t]. 
Since /(a) ^ 0, is monic — up to the scalar factor /(a) — of degree d in t. 

Set H for the gcd of F and G in Q(a;)[t]. Since F is monic in t, H belongs to Q[x, t], and it is easy 
to check that gcd(/, g) = H{x, 1) (up to a scalar factor). 

The procedure runs as follows: First we observe that if / = J2a fa{^ ~ o)", then 

^= E ( E Ux-aDt''-^ 

0<k<d \a\=k 

(and the same holds with g and G) so the homogeneous components of / and g centered at a turn 
out to be the coefficients of the monomial expansion of F and G with respect to t. Then, we apply 
the subresultant algorithm GCDl to compute the multiple Q e Q\x, t] and the superflous factor 
q in Qfx] of their gcd H in Q(a;)[i]. Finally, we apply Subroutine g to eliminate divisions in the 
expression gcd(/, g) := H{x,l) = Q{x,l)/q{x). 

Let us decide the size of the sets of points we have to take to insure that the algorithm has an error 
probability bounded by 1/4: 

We are going to choose randomly a point a G Q". This same point a will be used in each step we 
need a random point. 

The first condition the point a must satisfy so that the algorithm computes a greatest common 
divisor of / and g is that /(a) / 0. 
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Subroutine 3: Greatest Common Divisor 



procedure GCD{f,g,x,d) 

# fiQ & Q[a;i, • • • , Xn\ of degrees bounded by rf; 

^ X 1= {x\^ • • • ) ^n)': 

# The procedure returns h := gcd(/, 5). 

1. a := Random(n, 8 + 1)); 

2. if /(a) = then 

3. return( "error"); 

4. else 

5. (/o, ■■■,fd) ■■= Expand(/, x, a, d); 

6. (go,---,fi'd) :=Expand(fif,a;,a,d); 

7. e := 0; 

8. while geia) — and e < c? do 

9. e:=e + l; 

10. od; 

11. ife = d+lthen 

12. return(/); 

13. else 

14. := Eto fk t"""' and G := ELo 5fc t^"'; 

15. (g,Q) := GCDl(i^,G,d,e); 

16. h := PolynomialDivision(g(a;), Q{x, 1), d, q{a)); 

17. return(/i); 

18. fi; 

19. fi; 

end. 



Then we use the point a to compute the degree of G in t. Finally it is used in the subresultant 
algorithm to compute the degree of the gcd (by deciding whether certain determinants are zero or 
not). Checking the degree of G involves testing an r7-variate polynomial of degree bounded by d 
(the coefficients of G as a polynomial in Q[a;] [t]) while checking the degree of the gcd involves testing 
n-variate polynomials of degree bounded by 2d^. 

Thus, applying Schwartz bound for the set [0,^)", the conditional probability p of success verifies 

Therefore, taking £ := 8d{d+l) insures that the error probability is bounded by 1/4. 

Now let us compute the worst-case complexity of the machine: 

The cost of simulating the random choices here is 0{n logd). Computing the homogeneous com- 
ponents of / and g centered at a and checking the exact degree of G (that is finding the first 
non-zero coefficient of G with respect to t) can be done within complexity 0{d^{d + L)). In Al- 
gorithm GCDl, to compute the degree of the gcd involves computing at most d + 1 determinants 
of Sylvester-type matrices of size at most 2d x 2d, that is adds at most {d + 1) 0{d'^) operations. 
Once we know this degree, computing the gcd by means of an adjoint adds 0{d'^) steps. That 



19 



Subroutine 4: Power Series Division 



procedure PowerSeries(n, m, d, ipra, - ■■ , Vm+d, '4'm, • ■ • , ipm+d) 

^ n G N is the number of variables, 

# m G No is the order of the denominator (p G Q[[a;]], 

# d G N is the degree of the quotient h := ip/ip G Q[a;], 

# the ifiS and Vi's are the graded parts of the power series and i/i respectively. 

# The procedure returns q := ffyf^ h G Q[x]. 

3. (Po, ■■■:Pd)--^ Expand(P, t, 0, d); 

return(q); 
end. 



is, the complexity of computing Q{x, t) (and q{x) which is the non-vanishing determinant) is of 
order (D{d'^{d^ + L)) while L{Q{x,t),q{x)) — 0{(f{(f + LY) since the computation of the degree 
does not intervene in the length. Applying Subroutine g at q{a) which is different from (if 
not, subroutine GCDl in line 15 would have returned error) we obtain a final complexity of or- 
der 0{n \ogd + cP{d + L[Q{x,t),q{x)))) = 0{n log d + d^ (d^ + L)) . □ 

The following procedure (Subroutine ^ computes the quotient — provided it is a polynomial of 
bounded degree — of two multivariate power series from their graded components up to a certain 
bound. 

Let ip — Oq a;" G Q[[xi, . . . , be a power series. For i G No we denote by (pi :— 

S|a|=z '^ct -^^ G Q[xi, . . . , Xn] the j-graded component of (p. Also we denote by ordiy9 the order 
of (p, that is the least i such that (pi ^ 0. 

Proposition 1.9 Let p,tp G Q[[xi, . . . , Xn\] be power series such that h :— ^p/ip G Q[a;i, . . . , x„] . 
Assume we are given m := oidip, d > degh, and that the i-graded parts of p> and ip for i — 
m, . . . ,m + d are encoded by sip's of lengths bounded by L. 

Then PowerSeries (Subroutine ^ computes q := ^f^^ h within complexity 0{d^L). 
Proof - Set 

OC OO 

$(x, t) := p,{tx) = J2 i\ *(^: i^itx) = ^^(^) e ^ Q(a;)[M]. 

i=0 i=0 

Also set H := h{tx) G Q(a;)[t]. 

We first observe that ord$ = m, and thus ord^P >r7ias5'/$ = i7G Q(a;)[t] is a polynomial. 
Hence the following identity holds in Q(a;)[[t]]: 

H = - ^ — — (i- ^i^i^M!!^ ^ ^ V ~ 

Thus, for 



(J) pn in \ in I pn ^-^ 

^ t P-m \ Vm / t Pm 



\i=0 / \i=0 j=l 
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we have that tpf^^ H = P (mod Let P — Pi be the monomial expansion of P. Then 

V?^^ H = J2i=o ^ ^^'^ degree of H with respect to t is bomided by d. Hence (fif^^ h = J2i=o ^i- 
The stated complexity is obtained as follows: 

We compute a sip encoding of v := X^iLo?/'^ within complexity 0{d). We compute P as 
Yl'i=o'4'in+ii'' times v{Lpm, — Yl'j=i fm+jt^) withiu complexity 0{dL). We compute the expansion 
of P with respect to t up to degree d within complexity 0{d'^ L). Finally we compute q as Ylt^o Pi- 
The total complexity is of order 0{d^ L). □ 

Remark 1.10 In case that, in addition, we are given b £ Q" such that (pm{b) ^ 0, we can directly 
apply Subroutine^ to compute the quotient polynomial h within total complexity 0{d^L). 



2 The representation of the Chow form 



This section presents an algorithm for the computation of the Chow form of an equidimensional 
variety from a 0-dimensional fiber and a set of local equations at a neighborhood of this fiber. This is 
the key step in our general algorithm (see Section although it has independent interest: it shows 
that the Chow form and the geometric resolution are — up to a polynomial time computation — 



equivalent representations of a variety (see Subsection 3.1). As a further application, we give a 



non-trivial upper bound for the length of a sip representation of the Chow form (Corollary 2.12). 
In order to state the result, we need the following definitions: 

Definition 2.1 Let C P" be an equidimensional variety of dimension r. 

We say that fr+i, . . . , fn G ^ system of local equations at ^ lE V if the polynomials 

fr+i, . ■ ■ , fn generate I(V) at some neighborhood of i.e. I{V)(^ = (/r+i, . . . , /„)^ (where the 
subscript ^ denotes localization at the ideal of the point 

If Z is a subset of V , we say that /r+i, ■■•,/« £ liV) is a system of local equations (of V) at Z if 
it is a system of local equations at every ^ G Z. 

The existence of a system of local equations at a point ^ G V" implies that (C[a;]//(1/))^ is Cohen- 
Macaulay and thus, by Thm. 18.15], for /r+i, ••■,/« G I{V) to be a system of local equations 
at ^ is equivalent to the fact that the Jacobian matrix of this system has maximal rank n — r at ^. 

Definition 2.2 Let Z C A" be a 0-dimensional variety of cardinality D. A geometric resolution 
of Z consists of an affine linear form ^ = cq + ciXi + • • • + c„.x„ G Q[xi, . . . , a;„] and of polynomials 
p G Q[t] and v = (vi, . . . ,Vn) E such that: 

• The affine linear form £ is a primitive element of Z , that is £(^) ^ £{0 for all C ¥^ C i^i Z . 

• The polynomial p is monic of degree D and p{£{S,)) — for all ^ G Z ; that is, p is the minimal 
polynomial of £ over Z. 

• degVi < D — 1, I < i < n, and Z — {v{ri) ; rj E C / p{ri) = 0}; that is, v parametrizes Z by 
the zeroes of p. 

Observe that the minimal polynomial p and the parametrization v are uniquely determined by the 
variety Z and the afhnc linear form £. We say that (p, v) is the geometric resolution of Z associated 
to £. 

In case Z C P" is a zero-dimensional projective variety which satisfies that none of its points lie 
in the hyperplane {xq = 0}, Z can be identified to a 0-dimensional affine variety Z^^, the image 
of Z under the rational map P" A" defined by {xq : ■•• : x^) i— > {xi/xq, . . . ^x^/xq). By 
a geometric resolution of Z we then understand a geometric resolution — as defined before — of 
the affine variety Z^^ C A". In homogeneous coordinates, the definition of geometric resolution 
states that the homogeneized linear form satisfies £^/xo{£,) £'''/xo{£^') for all ^ 7^ ^' in Z. The 



21 



polynomial p is then the minimal monic polynomial of (i/xf) over Z^^ . On the other hand, v defines 
a parametrization V{p) ^ Z, r] ^-^ [\ : vi{'i]) : • • • ; v„(?7)). 

Now, we are able to state the lemma: 

Main Lemma 2.3 Let V C P" be an equidimensional variety of dimension r and degree D which 



satisfies Assumption l.i. Set Z :— V Cl V{xi, . . . ,Xr), and let p 6 Q[t] and v 6 Q[i]" be a 
given geometric resolution of Z associated to a given affine linear form i G Q[xi, . . . , Xn\- Let 
fr+i, ■ ■ ■ , fn S li^) be a system of local equations at Z. Assume that fr+i, . . . , fn have degrees 
bounded by d and are encoded by sip's of length bounded by L. 

Then there is a deterministic algorithm (Procedure ChowForm (Subroutine^ below) which computes 
Chv within complexity 0{r^ \og2{rD)n'^ d"^ D^^L). 



In Subsection 2.3 wc present the complete proof of the correctness of the algorithm and its complexity 



estimate. The algorithm is essentially based on a new Poisson type product formula for the Chow 



form (see Proposition 2.5 below), which describes the Chow form as a quotient of products of norms 
of certain polynomials. We interpret this expression as a quotient of two power series, which can be 
approximated with the aid of a symbolic version of Newton's algorithm. Finally we apply Procedure 
PowerSeries (Subroutine || above) to compute the Chow form from the approximate quotient. 



2.1 Newton's algorithm 

In this subsection we present a symbolic version of Newton's algorithm for the approximation of 
roots of equations. Newton's algorithm is a widely used tool for polynomial equation solving. The 
situation in the present work is not much different from that in e.g. p7| . Hence we just describe 
this procedure in order to adapt it to our setting and to estimate its complexity; its correctness follows 
directly from Section 2] and the arguments therein. 

First, we state the situation in which Newton's algorithm is applied. 

Let C A"* X A" be an equidimensional variety of dimension r such that the projection map 
TT : W ^ A"^ is dominant, that is, the image tt(W) is a Zariski dense set. 

Set A := Q[ti, ...,tr] = Q[A''] and let K be its fraction field. Also let B := Q[W] and set L := 
K ®A B. Then L is a finite if-algebra, and its dimension D := [L : K] — that is the degree of tt — 
equals the maximum cardinality of the finite fibers of n (p9|. Prop. 1]). 

The norm N^(/i) 6 if of a polynomial h £ A[xi, . . . , Xn] is defined as the determinant of the 
iiT-linear map L ^ L defined hy f ^ h f . 

Let I{Wy denote the extension of the ideal L{W) to the polynomial ring K[xi, . . . ,Xn], and set 
:= V{I{Wy) C A'^{K), which is a 0-dimensional variety of degree D. Then 

^Ah)= n '^w- 

We also denote this norm by Nw^{h) when the projection map is clear from the context. 

In different steps, we will be given a polynomial h and an equidimensional variety W and our aim 
will be to compute an approximation of Nvyc (/i) . The input of the procedure will be the polynomial 
h, a geometric resolution of a 0-dimensional fiber of tt and local equations at this fiber. 

Let Fi, . . . , F„ e I{W) C A[xi, ...,Xn\ and set F := {Fi, . . . , F„). Let 

'^^-(tt) €^[xi,...,a:„]"><" 
\oxj J .,■<„ 

be the Jacobian matrix of F with respect to the dependent variables xi, . . . ,a;„, and let :— 
\Jf\ G A\x\^ . . . , Xn] be its Jacobian determinant. 

Let Z C A" be such that 7r~^(0) = {0} x Z. We assume that Z is a 0-dimensional variety of 
cardinality D and that Jp is non-sing ular at 7r-i(0), that is, Af(0,O 7^ for all ^ G Z. Observe 
that this means that i^(0, x) is a system of local equations at Z. 
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Under our assumptions, the elements in can also be considered as power series: For ^ ^ Z, the 
fact that Ap{Q, ^) 7^ implies that there exists a unique 7^ G C[[ti, . . . , tr]]" such that: 

7^(0) =e and F{tu...,tr,li)=0. 

(See, for example, ||, Ch. 3, Section 4.5, Cor. 2].) It follows that /(^i, . . . , t^, 7^) = for all 
/ G I{W) as F is a system of local equations at and so this also holds for all / e I{W'^) — I{WY. 
Hence 7^ G and, as #Z = = D, we conclude that the correspondence 

Z^W , C ^ 7? 

is one-to-one. In particular, since (h) is the determinant of a matrix in Q(ii, . . . ,tr), Nw {h) G 
C[[ti, . . . , U]] n Q(ti, ...,U)C Q[[ti, . . . , i„]]. 

The given data — the description of the fiber and its local equations — suffices to determine 
uniquely, and in particular allows to compute a rational function q which approximates the norm 
Nvye (h) to any given precision k (we understand by this that both Taylor expansions coincide up to 
degree k, that is Nvk<=(/i) = q mod {ti, . . . ,trY'^^). The rational function q can be obtained by a 
procedure based on an iterative application of the Newton operator. This operator, defined as 

:= x' ^ Jf(x)-^ F(xY ^K{xT''-\ 

enables us to approximate the points in W from the points in the fiber Z. If we set TVj,™-' G K{xY^^ 
for the TO-times iteration of A/i?, then, for every f G 

a4")(0^7« mod (ii,...,t,)2" 

(see Section 2]). 

Procedure NumDenNcwton (Subroutine ^ computes polynomials 5^™^ , . . . , gi" 
such that 

AA(") = (5^V/o^"\...,5.^V/o^"^ 

Herein, Homog(/, d) is a procedure which computes the homogeneization of the polynomial / up to 
degree d > deg(/), JacobianMatrix(i^, x) is a procedure which constructs the Jacobian matrix with 
respect to the variables x associated to the system of polynomials F and Adjoint(Af ) is a procedure 
which computes the adjoint of the matrix M. For the correctness and complexity of the whole 
procedure, see Lemma 30]. 
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Subroutine 5: Computation of numerators and denominators for the Newton operator 



procedure NumDcnNewton(i^, rt, x, d, m) 

# Fi,. . . ,Fn e A[xi, . . . ,Xn] such that Jf{x) ^ 0; 

# n is the number of dependent variables a;, 

# d is an upper bound for the degrees of the polynomials Fi, . . . , i^„, 

# m is the number of iterations to be computed. 

# The procedure returns g["'\ gt\ ft^ such that Mp{x)("^') = (ff^V/o^"^ • • ■ , gt^/ft^) 

1- J^F ■= JacobianMatrix(F, x) ; 

2. Af := |(Jf)|; 

3. A := Adjoint(J>); 

4. z/ := + 1; 

5. for i from 1 to n do 

6. ,gf ^ := Afi-, - E"=i ^^i/jS 

7. G, :=Homog(.g|i\i.); 

8. od; 

9. f!,'^ := Af; 

10. Fo := Homog(AF,z^); 

11. for k from 2 to m do 

12. for i from 1 to n do 

14. od; 

lO- Jo ■— -f^oUo '5l ,---i9n ), 

16. od; 

return(5j™\ . . .,gi"'\ /,$™^); 
end. 



We summarize the procedure that approximates the norm of a given polynomial h in Procedure 
Norm (Subroutine ^). Herein, CompanionMatrix is the procedure which constructs the companion 
matrix of a given univariate polynomial. We keep our notations. 

Lemma 2.4 Let notations be as before. Assume that h, Fi, . . . , F^ £ A[xi, . . . , a;„] are polynomials 
encoded by sip's of length bounded by L such that degh < S and deg(-Fi) < d, 1 < i < n. 
Then Norm (Subroutine ^ computes f,g E A with /(O) ^ such that g/ f approximates Nn/c(/i) 
with precision k, within complexity ©((logj k) rJ S^cP' D'^L). 

Proof - For the correctness of the algorithm we refer to Section 2] and the arguments given 
there. Now, we estimate its complexity: 

First, the complexity of Subroutine ^ applied to our situation is of order C'((log2 n) n'^d^L) (see pO| , 
Lemma 30]). 

Then, the algorithm computes the matrices Vj{Mp) (1 < j < n) with complexity of order 0{nD^) 
(note that, as the companion matrix is very sparse, the multiplication by Mp can be done with 
complexity 0{D^)). Now, the matrices gi{v{Mp)) {1 < i < n) and Mq := fo{v{Mp)) 
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Subroutine 6: Approximation of the Norm 



procedure Norm(ft,, S,n, x,p,v, F, d, k) 

=ff h Cz A[xi, . . . , Xn] is the polynomial whose norm we want to approximate, 

# 5 is an upper bound for the degree of h, 

# n G N is the number of dependent variables x, 

# p e Q[t], V e Q[t]" is a given geometric resolution of Z, 

# F = (i^i, . . . , Fn) is a vector of polynomials in I(W) such that Ai?(0, ^) 7^ for all ^ € Z, 

# d is an upper bound for the degrees of the polynomials Fi, . . . ,Fn, 

# K e N is the desired level of precision. 

^ The procedure returns f,g & A with /(O) 7^ such that g/ f approximates the norm 'Hwih) 

# with precision k. 

1. m [log2(K+ 1)1; 

2. ((71, . . . , (7„, /o) := NumDenNewton(F, n, d, m); 

3. Mp := CompanionMatrix(p); 

4. for i from 1 to n do 

5. M, := g^{v{Mp))■ 

6. od; 

7. Mo /o(«(Mp)); 

8. H := Homog(/i,^); 

9. M ■.^H{Mo,Mi,...,Mn); 

10. /:= |Mo|^; 

11. (7:= \M\; 

return(/,g); 
end. 



are obtained within complexity 0{{\og2K)'n7(fD^L). As /i is encoded by a sip of length L, its 
homogeneous components up to degree 6 are encoded by sip's of length 0{S^L). Therefore, the 
complexity of the computation of M is of order 0{6^LD^ + (log2 k) 11^ d?D^L). 
Finally, / and g can be computed within complexity 0{D^ + (log2 6) + (log2 n) rJcPD^L) and 
0{D'^ + + (log2 k) v7(fD'^L) respectively. □ 

2.2 A product formula 

In what follows, we establish a product formula for the Chow form of an afhne variety. This formula 
is an analogue of the classical Poisson formula for the resultant (|l^, Ch. 3, Thm. 3.4]). It describes, 
under certain conditions, the Chow form in a recursive manner. 

Let V C A" be an equidimensional affine variety of dimension r and degree D which satisfies 



Assumption 1.2. Let Uq, . . . ,Ur he r + 1 groups of n + 1 variables each and let Li := Uio + Un xi + 
••• + UinXn, < i < r, be the affine linear forms associated to these groups of variables. Set 
K := Q{Uo, . . . , Ur^i) and let I{Vy denote the extension of the ideal of V to the polynomial ring 
K[xi, . . . , Xn] (or to any other ring extension of Q[a;i, . . . , Xn] which will be clear from the context). 
We also set V'^ :— V{I{VY)r\V{Lo, . . . , Lr-i) C A" (if ), which is a 0-dimcnsional variety of degree 
D. 
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For < i < r, we set Vi := y n V{xi+i, . . . ,Xr) C A", which is an equidimensional variety 
of dimension i and degree D as V satisfies Assumption 1.2. Observe that these varieties satisfy 
Assumption 1.2 as weh. 
Let Ki Q{Uo, . . . , C/j-i) ^ K and set 

V{I{V,y ) n 1/(^0, . . . , L,_i) C K^(Ki). 

Observe that is also a 0-dimensional variety of degree D. 

Under these notations we have that V^'' = Vb, Kr ~ K and = . 

Proposition 2.5 (Product formula) Let V C A" he an equidimensional variety of dimension r 
which satisfies Assumption Let notations be as in the previous paragraph. Then 



Y{Chyo{Ui) 

ChviUo, ...,[/.)= '-^ ^ e Q{Uo, Ur^i)[Ur]. 

n Chyo(ei) 
1=1 

The proof of this fact is based on the foUowing lemma: 

Lemma 2.6 Let V C A" be an equidimensional variety of dimension r. Let Ty G Q[C/o, ■ • ■ , Vr\ 
and J-'yo e ^^[J/r] be Chow forms of V and V'^ respectively. Then there exists X G K* such that 

Ty — A Tya . 

Proof.- As before, we denote by l^YY the extension of the ideal L{y") to a ring extension of 
Q[a:;o, . . . ,a;„] which will be clear from the context. Let J7'™, < i < r, denote the group of n 
variables Ui \ {Uio). We consider the map 

Q[{/o, . . . , c/,][xi, . . . , x^\i(i{vr + (Lo, . . . , U)) ^ qK'", . . . , c/M[xi, . . . , x^ynyr 

defined by f/io — (CA i 2:1 + ■ • • + a;„), C/y 1— > Uij and for < i < r, 1 < J < 

As it is a ring isomorphism, I{VY + (^O: ■ • • j ^r) is a radical ideal. Following notations in Subsection 



LI, it follows that this ideal coincides with the defining ideal of the incidence variety <^y and, 



therefore, 

{Ty) = {I{Vr + {Lo,...,Lr)) n Q[{/o,...,C/.]. 

Similarly 

{Ty,) - + (L,)) n K[Ur]. 

We have that + {Lo,...,Lr) C /(V^")" + [Lr) and so (J^y) C [Ty^), that is, there exists 

A e if[i7r] \ {0} such that Ty = XTyo. As deg^^ Ty = deg V = degTyo, A is an element in K* . □ 



Proof of Proposition 2.5.- Let 1 < « < r. From Lemma |2.q , there exists \i £ K* such that 

Chy{Uo,...,U,)^X,Chyo{Ui). (1) 

Hence C/iy^ (fJoj ■ ■ ■ , C/j-i, e^) = A^ Chyo (ej. 

Now, it is easy to see that Chy^ -^ (Uq, . • . , Ui-i) divides Chy. {Uq, . . . , Ui-i, e;). From Assumption 



1.2, it follows that degl^_i = degT^ = D and, therefore, both polynomials have the same degree. 



Moreover, the normalization imposed to both Chow forms implies that they coincide. So 

Chy._^{UQ,...,Ui^i) = XiChyo(ei). (2) 
From Identities (|I]) and (||) we deduce that 

Chy{Uo, . . . ,Ui) Chyo{U^) 



Chy.^^{Uo, . . . ,Ui_i) Chyo{ei) 



(3) 
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Multiplying these identities for z = 1, . . . , r we obtain 



Chv{Uo,...,Ur) ^yr Chy^jUp, . . . ,Ui) ^^ ChypjUj) 

ChvoiUo) Chv^_,{Uo, C/.-i) Chvo{ei) 



which gives the formula stated in Proposition 2.5. □ 



Observe that Chyo{Ui) = Yl^izyo ^lij) = Nyp{Li) and Chya{ei) = 1176^0 ^*('^) ^ Ny.o(a;i), which 
implies that the Chow form of V can also be presented as the quotient of two products of norms of 
polynomials: 

Corollary 2.7 Let notations and assumptions he as before. Then 

Chv{Uo. ...,Ur) = ^ e QiUo, Ur-l)[Ur]. 

n Nyo(x,) 

This formula enables us to compute Chy as a quotient of power series. To do so, we are going to 
prove a technical lemma first. 

Lemma 2.8 Let V C A" be an equidimensional variety of dimension r which satisfies Assumption 
Assume that V is Cohen- Macaulay at every point of Z := V H V{xi^ . . . , Xr)- Then, the ideal 
I{V) + {xi, . . . ,Xr) C C[a;i, . . . ,a;„] is radical. 

Proof - Let V denote the projective closure of V C A" ^ P". Let Z := V n V{xi, . . .,Xr). The 
fact that #Z = deg V imphes that Z = V Ci V{xi, . . . , Xr). 

Take ^ G Z and let be the primary component of the ideal I{V) + {xi, . . . , Xr) C C[a;i, . . . , Xn] 
which corresponds to ^. 

We consider the length £(V, V{xi, . . . , Xr); which under our assumptions can be defined as 

£(F,F(xi,...,x,);0 =dimcC[xi,...,a;„]/Q5. (4) 
By a suitable version of Bezout theorem (see |5^, Prop. 3.30]) 

Y,HV,V{xi,...,Xr);0<<iegV. 

On the other hand, as i{V , V{xi, . . . , Xr)', is a positive integer for each £, G Z and #Z — deg^, 
it follows that 

Then liV , V{xi, . . . , Xr); = 1 fo^' ^ ^ and so (Q) implies that = {xi — ^i, . . . , a;„ — 
which is a prime ideal. 

As / := I{V) + {xi, . . . , Xr) is 0-dimensional it has no embedded components. Hence / = is a 

radical ideal. □ 

The following corollary shows that the coordinates of all the points in belong to the subring 
C[[[/o — ei, . . . , Ui-i — Ci]] n Ki, and that there is a one-to-one correspondence between the points 
oi Z V n V{xi, . . . , Xr) and the points of V^. 



Corollary 2.9 Let notations and assumptions be as in Lemma 2.1 and before. Let < i < r and 

^ E Z. Then there exists a unique 7^*' G C[[C/o ~ ei,...,C/i_i — ej]" such that 7_|'^ £ and 
7^'^(ei,...,e,) = 
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Proof.- Suppose I(V) is generated by the polynomials hi, . . . , hf. Since we are in the conditions 
of the previous lemma, the Jacobian criterion jl^, Thm. 18.15] implies that the Jacobian matrix 
associated to the generators hi, . . . , ht, xi, . . . , Xr of the ideal I{V) + {xi, . . . , Xr) has maximal rank 
n at ^. In other words, there are n polynomials gi, . . . ,gn among hi, . . . , ht, xi, . . . , Xr such that the 
associated Jacobian determinant is non-zero. Now, as the rank of the Jacobian matrix of hi, ht at 
^ is bounded by the codimension n — r of V at we can assume w.l.o.g. that gi := Xi, . . . , gr ■= Xr- 
Let 

'% 

. dx. 



A := 



3 / r-]-l<.i.'j<.n 



Then A(0 7^ 



be the Jacobian determinant of gr+ii ■ ■ ■ , gn with respect to the variables Xr+i, ■ ■ ■ ■ 
since A coincides with the Jacobian determinant of the system gi, . . . , g„. 
On the other hand, let A^ £ Ki[xi, . . . , Xn] denote the Jacobian determinant of the system 
Lq,. . . ,Li_i,Xi+i, . . . ,Xr,gr+i, ■ ■ ■ ,gn- An casy verification shows that Ai{ei, . . . ,ei_i)(^) = A(^) 7^ 
0. The statement follows from the arguments in Subsection 2.1. □ 



Now, set 

r r 
^ ■.= l[fiv«{L^) eK[Ur], $ :=[|Nv;o(x,) & K* , 



i=0 i=l 



SO that, by Corollary |2.7], Chy := 

From Corollary U, f e Q[[Uo - ei, . . . , Ur-i - er]][Ur] and $ e Q[[Uo - ei, . . . , Ur-i - e^]]. 

The following lemma gives the order of the denominator ^ at E :— {ei, . . . , Cr) & A*" together 
with its graded component of lowest degree: 

Lemma 2.10 Let notations be as in the previous paragraph and let D :— deg V^. Then ord _b($) — 
rD and its graded component of degree rD is 

r 

^rD = ±l[Chv„{Ur-l). 

1=1 

Proof - Clearly, ordB($) = ELi ord£(Nv^^o(a;,)). 

Let 1 < i < n. Recall that Nyo(a;i) = C/iyo(ei). From Identity (|^) in the proof of Proposition |2.5| 
we have: 

Chyolei) ChvAUo, ■ ■ ■ -.Ui) ^ Chyo{Ui) Chvi_i{Uo, . . . , Ui-i). 

As C/iy.o(eo) = 1, then Chyo{ei) Chvi{Uo, ■ ■ ■ , J7i_i, cq) = Chvi_^{Uo, ■ ■ ■ , Ui-i). 
We also have that 

Chv,{ei, . . . ,ei,eo) = ±C/iy, (ep, ei, . . . ,6^) = ±1. 

This shows that Chvi {Uq, . . . , Ui^i, cq) is invertible in Q[[[/o — ei, . . . , Ui-i — ei]] and, therefore, if 
TO := ord£;(C/iyo(ei)), 



C/iyo(ei) = ±Chvi_j^{Ua, . . . ,Ui-i) mod {Uq - ei, . . . ,Ui-i - ei) 



m+l 



By Lemma 2.6, there exists Ai_i G Q{Uo, . . . , Ui-2) \ {0} such that 

Chvi_i(Uo, . . . ,Ui^i) = \i-iChvo_^{Ui-i). 

As ChyQ ^{Ui-i) is a homogeneous polynomial of degree D in the group of variables Ui^i and does 
not depend on Ui-n, it is also homogeneous as a polynomial expanded in Ui-i — ei. Then, the 
order of Chvi_i at with respect to the group of variables Ui-i equals D. On the other hand, we 
have that Chvi_i{ei, ■ ■ ■ ,ei-i,Ui-i) = ±Chvo{Ui-i) ^ 0. This implies that the series Chvi_i in 
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/o — ei, . . . , Ui-i — Ci]] has a term of degree D depending only on the group of variables Ui-i — e^. 
We conclude that m = D and 

{'Nyo{Xi))D = {Chya{ej))D ±(C/ly._i([/o, . . . , U,^i))d = ±C/ly„(C/^_i). 

Therefore, ord_E(<i') = X]I=i 0''d-E(Nyo(xi)) = rD and the graded part of lowest degree of $ is 
2.3 The algorithm 

Here, we are going to put the previous results together in order to obtain the algorithm underlying 
Main Lemma |2.3| and to estimate its complexity. 



Let notations be as in Main Lemma 2.3. As we have already noted, the imposed conditions imply 
that both V and Z have no component in the hyperplane {xq = 0}. Hence V equals the projective 
closure of its afhne part Vxq := V \ {xq = 0} and so their both Chow forms coincide. Hence we 
concentrate w.l.o.g. on the afRne case. We use afSne coordinates and keep the notation of the 
previous subsection. 

From Corollary 2.7, we have that 

flNv.o(L,) 
Chv^'-^ . 



n Nyo(a;,) 

i=l 



Now, we approximate the norms appearing in this formula. 
Set 

V, V{I{V,r) n y(Lo, . . . , i^-i)) C A'("+i) X A". 

The map tt^ : Vi — > A'*^"+^) defined by {U,x) i— > C/ is dominant of degree D := degT^. We set 
Z ■= Vo = VnV{xi,...,Xr) C A" and let E, := (ei,.. .,6^) £ A*("+i). Then 

-.^ 7:r\E,) = {E,} X Z 

and so this fiber is a 0-dimensional variety of cardinality D. Furthermore, it is easy to check that 

Lo, . . . , ij-l, Xi+i, . . . ,Xr, fr+1, ■ ■ ■ ,fn G Q[Uo, ■ ■ Ui-i][xi, . . . , Xn] 

is a system of local equations of Vi at Zi. 

Since by definition, 'NT^.{xi) and N7r;(ii) coincide with Nyo(a;i) G Q[[C/o — ei, . . . , Ui-i — €{]] and 
Nyp(Li) G Q[[C/o — ei, . . . , J7i — Bi+i]] respectively, we can compute any approximation of the lat- 
ter applying Procedure Norm (Subroutine ||) modulo a change of variables (C/q, . . . , Ut) t-^ {Uq + 
ei, . . . ,Ui + Gi+i) (in order to center the series at 0). 

We multiply the computed approximations for < i < r to obtain rational functions ■(/' and Lp which 
approximate the power series 

r r 

■.= l[Nyo{L,) eQ[[Uo-ei,...,Ur-i~er]][Ur], <!>■.= l[l^yo{x,) eQ[[Uo-ei,...,Ur-i-er]]* 

1=0 i=l 

respectively. 

From these approximations, we compute the graded parts of $ and ^' of degrees between rD and 
{2r + l)D centered at (E, 0) e A(''+i)("+i) (where E :^ {a, . . . ,er) € A''("+i)) by applying Proce- 
dure GradedParts (see Subsection 1.4). 

By Lemma 2.10 , we have that ord (b_o)('J') = ord b($) = rD. We also have degChy = (r+l) D. We 
use this information together with the obtained graded parts in order to apply Procedure PowerSeries 
(Subroutine This yields a polynomial Q £ Q[C/oi ■ • • j Ur] such that 

Q^^i;+^)^+^Chy. 
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Subroutine 7: Chow form from a fiber and local equations 



procedure ChowForm(n, x,r, D,p,v, f,d) 

^ n is the number of variables x := {xi, . . . , Xn), 

# D are the dimension and the degree of V respectively, 

# ]3 G Q[t], V G Q[t]" is a given geometric resolution of the fiber Z, 

# / = (/r+i, ■ • • , /n) G Q[a;i, ■ • • , Xn]""'' is a system of local equations oi V at Z of degrees 

# bounded by d. 

# The procedure returns the normalized Chow form Chv- 
1. for i from 1 to r do 

2- ■= Norm(a::j, l,n,x,p,v, Lq, . . . ,Li„i,Xi+i, . . .,Xr,fr+i, . . . , /„, d, (2 r + 1) D); 

3. od; 

4. for i from to r do 

5. (V'f \ V'i^'') Norm(Li, l,n,x,p,v, Lq, . . .,Li^i,Xi+i, ...,Xr, /r+i, ...,/„, d, (2 r + 1) £)); 

6. od; 

9. ($0, • • • , *(2r+i) d) GradedParts((^(i), </7(2), (ei, . . . , e,), (2 r + 1) D); 

10. (*o,...,*(2r+i)D) := GradcdParts(V'W,i^(2)^ (ei, . . . , e,, 0), (2 r + 1) i?); 

11. Q := PowerSeries((r + 1) (n + 1), r £>, (r + 1) D, . . . , «'(2,.+i) d, *rr>, • • ■ , *(2r+i) zj); 

12. C/iv := PolynomialDivision(Q, (r + l)D, (eo, . . . , eo)); 
return(C/iv); 

end. 



Again, from Lemma ^.IC , the denominator <i>^'^^''^^^ does not vanish at (eo,...,eo) G 

p^r(n+i) -y^g apply Proccdurc PolynomialDivision (Subroutine ||) to the polynomials Q and 
and the point . 

We summarize this procedure in Procedure ChowForm (Subroutine^) which computes the Chow 



form of an afSne equidimensional variety V satisfying Assumption 1.2 



Proof of Main Lemma 2.c As we have already observed, we may suppose w.l.o.g that V is an 
affine variety and that the polynomials /r+i, ■ ■ ■ , fn are in Q[xi, . . . , a;„]. 

We apply Procedure ChowForm (Subroutine |^) to V in order to compute its normalized Chow form. 
The correctness of this procedure follows from our previous analysis. The announced complexity is 
a consequence of the complexity of the subroutines we call during this procedure: 

• By Lemma |2.4| , the complexity of lines |^ to ^ is of order 0{r \0g2ir D) in7 L) . The 
products in lines ^ and || do not change this estimate. 

• The computation of the graded parts in lines p| and has complexity 0{r^ log2(?' D) d? D^L). 

• Finally, the subroutines PowerSeries and PolynomialDivision in lines |ll| and [l^ add complexity 
0{r^ \og2{rD)rJd'^D^^L). 

We conclude that the overall complexity is 0{r^ ^og2{r D) n7 d^ D^^L). □ 
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Remark 2.11 The Chow form Chy is the numerator o/Nyo(Lr) by Lemma 2.6. Unfortunately 
this norm is a rational function, due to the fact that the map tt^ is not finite but just dominant. 
The product formula is the tool which allows to overcome this difficulty, as it gives an expression for 
Chv without any extraneous denominator. 



We directly derive the following estimate for the length of a sip representation of the Chow form of 
an equidimensional variety: 

Corollary 2.12 Let V C P" be an equidimensional variety of dimension r and degree D. Let 
fr+i, ■ . ■ , fn G ^(^) be a system of local equations at a dense open subset of V , encoded by sip's of 
length bounded by L. Then, if d := max{deg(/i) : r + 1 < i < n}, we have 

L{Fv) < 0{r^ \og^{rD)n''d^D^^L). 

Proof.- Let ti, . . . ,lr G Q[a^O: ■ ■ • i Xn] be linear forms such that Z -.^ V Ci V{£i, . . . ,£r) is a 0- 
dimensional variety of cardinality D. We can choose these linear forms so that Z lies in the dense 
open subset where /r+i, ■ . . , fn is a system of local equations. 

Furthermore, let £o, £r+i, ■ ■ ■ ,£n be linear forms which complete the previous ones to a change of 
variables such that Z n {^o = 0} = 0- 

Then V satisfies Assumption |l.2| with respect to these variables, and the statement follows directly 



from Main Lemma 2.3. □ 



3 The computation of the Chow form 



We devote this section to the description and complexity analysis of the algorithm underlying The- 
orem 1. The first subsections gather some results which lead to the proof of the theorem. 



3.1 Geometric resolutions 



Geometric resolutions where first introduced in the works of Kronecker and Konig in the last years 
of the XlXth century. Nowadays they are widely used in computer algebra, especially in the 0- 
dimensional case, but there are also important applications in the general case. We refer to |^ for 
a complete historical account. 

In what follows we show that we can compute any — sufficiently generic — geometric resolution of an 
equidimensional variety from a Chow form in polynomial time. This computation and the procedure 
described in Section || imply that, from the point of view of complexity. Chow forms and geometric 
resolutions are equivalent representations of an equidimensional variety. 

Let V C A" be an equidimensional affine variety of dimension r and degree D. For < i < r, let Li 
denote as usual the generic affine forms. Let Ci G Q"+^. We set 

4 := Li{ci) = Cm + Cii xi ^ h x„ e Q[xi, . . .,x„]. 

We assume that the projection map 

T^{ei,...,er) '■'^ > X ^ {£i{x), . . . Jrix)) 

is finite, that is the affine linear forms £i, . . . ,£r are in Noether position with respect to V. Let 
yi, . . . ,yr be new variables. Set 

K ■.= Q{yi,...,yr) , L = Q{£i, . . . ,£r) (»q[e,_,e,^] Q[V] 

and consider the morphism 

K ^ L , £i. 
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Then K ^ L is a finite extension of degree [L : K] < D. We assume furtliermore tliat £o is a 
primitive element of this extension, that is L = -R'[^o]- 

Then the geometric resolution of V associated to I := {£o, . . . , £r) is the pair 

p -.^ pvj e K[t] , w := wvj e K[t]" 

where p is the monic minimal polynomial of £o with respect to the extension K ^ L, and w = 
{wi, . . . , Wn) verifies degwi < [L : K] and p'{£o) Xi = Wi{£Q) e L for 1 < i < n, where p' := dp/dt. 
These polynomials are uniquely determined and because of the Noether position assumption, we 
have that p,Wi lie in fact in Q[yi, . . . ,yr][t], see e.g. j27[ Section 3.2]. 

A geometric resolution gives a parametrization of a dense open set of V in terms of the points of a 
hypersurface in A'"+^: there is a map 

V{p{t, yi, . . . , yr)) \ V{p\t, yi, . . . , yr)) ^ V \ Vip'{£o{x), £i{x),..., £r{x))) 

(t,yi,...,yrj ^ —{t,yi,...,yr}. 
P 

Note that, in case the considered variety is 0-dimensional, this definition of geometric resolution 
essentially coincides with the one given in Section |[ the passage from one to the other can be made 
by considering the resultant with respect to the variable t between p and p' . 

The following construction shows that the geometric resolution associated to the generic affine linear 
forms Lq, . . . , Lr can be expressed in terms of the characteristic polynomial of the variety, and hence 
in terms of the Chow form: 

Let Uq, . . . ,Ur he r + 1 groups of n + 1 variables which correspond to the coordinate functions of 
j^(r+i)(n+i) g^j-^^ jg-j. rp j-j,^^ ^ J,^^J ^ group of r + 1 variables which correspond to the coordinate 

functions of A''+^. We recall that a characteristic polynomial Vv £ Q[C^i • ■ • j ^^r][7o, . . . , TV] of V 
is defined as any defining equation of the Zariski closure of the image of the map 

: A('-+i)("+i) xV-^ A('^+i)("+i' X A-'+i, (uq, . . . , u,; ^ K, Lo{uo, ■ ■ • , ^^K, 0) 

which is a hypersurface. This is a multihomogeneous polynomial of degree D in each group of 
variables Ui U {Ti}. Its degree in the group of variables T is also bounded by D. 

A characteristic polynomial of V can be derived from a Chow form For 1 < i < r we set 

Ci := {Uio - Ti, Uii, Uin). Then 

= (-l)^J^v(Co,...,Cr) (5) 

is a characteristic polynomial of V . 

Set Vv '■= <idTq+- • -+00 for the expansion of Vv with respect to Tq. Then od lies in Q[f/i, . . . , [/r]\ 
{0}, and in fact it coincides with the coefficient of f/J^ in Tv, that is 

aoiUi, ...,Ur)^ ^y(eo, Ui,..., Ur). 

In case V satisfies Assumption |l.2| , we define the characteristic polynomial of V as 

(-l)^C/lv(Co,...,Cr-) 

where Chv is the normalized Chow form of V . We refer to Section 2.3.1] for further details as 
well as for the proof of the stated facts. 

Lemma 3.1 Let V C A" he an equidimensional variety of dimension r and degree D. Let Uq, . . . ,Ur 
be r+1 groups ofn+l variables and let Lq, . . . , Lr be the generic affine forms associated to Uq, . ■ ■ , Ur- 
Set E :— Q(C/o, • . • , Ur) and let denote the Zariski closure of V in A"(i?). Let To, . . . ,Tr be new 
indeterminates. 



32 



Then the geometric resolution of associated to Lq, ■ ■ ■ , Lr is given by 



P:^^eE[T,,...,Tr][To] , W := -— ( . . . , e E[T,, . . . ,Tr][ToY^ 

ao an \oUoi oUonJ 

where Vy is a characteristic polynomial of V and ajj is the leading coefficient of Vy with respect to 
To. 

Proof.- Using the fact that the extended ideal C E[xi, . . . , Xn] is radical, it is easy to check 

that I{V°) — I{Vy. Consider then the morphism 

A:=E[Ti,...,Tr]^B:=E[xi,...,x„]/IiVr , T,^L,{U„x). 

Our first aim is to prove that this is an integral inclusion, or in other words, that the projection 
map TT(^Li,....Lr) ■ ^ A'^{E) is finite. 
By definition 

P{Uo,...,Ur)iLoiUo,x),...,Lr{Ur,x))=0 mod /(F) ^((7) [a;] . (6) 

Specializing Uo by the {i + l)-th element of the canonical basis Ci in this identity, we deduce that 
P{ei,Ui, . . . ,Ur){To,Ti, . . . ,Tr) G ^[Tq] is a monic equation for Xi for i = Therefore 
^ ^ S is an integral extension. 
Set /C := E{Ti, . . . ,Tr) and C:=JC®a B. 

It is immediate that P :— Vv Io,d is a monic polynomial equation for Lq with respect to the extension 
K. ^ L. As ^ ^ S is an integral extension, from the definition of Vy we deduce that P is the 
minimal monic polynomial of Lq. This implies that : /C] = Z? and that Lq is a primitive element 
of this extension. 

Write 

Q(f7, x) Vv{Uo, . . . , Ur){Lo{Uo, x),..., L,(C/„ x)) = Y.bp 

P 

with 6/3 e Q[C/i , . . . , [/r] , . . . , a;„] 



As bp e I{VY C Q[C/i,...,C/.][xi,...,a;„] for aU /3, G I{Vy C Q[C/o, . . . , i7.][xi, . . . , 

dQ(U,x) 

for i = 1, . . . , n. Therefore, — - — ^ — = in £ for 1 < i < n. Then the chain rule implies that the 

oUoi 

identity 

r)P f)P 

-f{u, m x)) = -w^iu, m ^)) 

Ol OUoi 

holds in C and the lemma follows. □ 

Now we show how a particular geometric resolution can be obtained by direct specialization of the 
generic one. 

Using the same notation as in the beginning of this subsection, we will assume that Vr\V{ii, . . . ,£r) 
is a 0-dimensional variety of cardinality D. This condition is satisfied provided that ii, . . . ,£r are 
generic enough [Bq, Prop. 4.5]. After a linear change of variables, we may assume w.l.o.g. that 



£i = Xi for i = 1, . . . , r, so that the stated condition is Assumption 1.2 
Thus, for the rest of this section we fix the following notations: 

Z ■.^Vr]V{xi,...,Xr) , K ■.= Q{xi,...,Xr) , L '.^ K i 

We also assume that £q = Lq{co, x) G Q[xi, . . . , Xn] separates the points of Z. This is also a generic 
condition: if we set p := discrj-Q Vz G Q[t^o] \ {0}, this condition is satisfied provided that p{co) ^ 0. 
These two conditions ensure the existence of the associated geometric resolution of V: 
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Lemma 3.2 Let V C A" be an equidimensional variety of dimension r and degree D which satisfies 



Assumption l.i. Let Ho :— Lq{cq,x) € Q[xi, . . . , a;„] be an affine linear form which separates the 
points of Z . 

Then the projection map tt : V A'' , tt{x) — (xi, . . . , Xr) is finite and £o is a primitive element of 
the extension K ^ L. The geometric resolution of V associated to £ := {£o,Xi, . . . ,Xr) is given by 

p := Vvico,ei, . . . ,er)it,xi, . . . ,Xr) G . . . , x,.][t], 

w ■■= - f (co,ei,...,er)(t,a;i,...,a:r-) G Q[xi, . . . , a;^] [t]", 

where Vv is the normalized characteristic polynomial of V . 

Proof.- The fact that tt is finite foUows from Lemma 2.f4]. On the other hand, the normalization 
imposed on Vv implies that 

PoW :=p(i,0,...,0) =7'v(co,ei,...,e,)(t,0,...,0) £ Q[t] 

is a monic — and thus non-zero — polynomial of degree D which vanishes on i?o(0 for all ^ G Z . The 
hypothesis that £{) separates the points of Z implies that po is the minimal polynomial of over Z; 
in particular it is a squarefree polynomial of degree D and so, as p is monic, 

7^ discrpo — (discr(p)(0, ... ,0). 

In particular, discr^p ^ and thus p is also a squarefree polynomial which annihilates £q over V . 
Now, as the map tt is finite, the minimal polynomial m^,-, G K[t] of lies in Q[a;i, . . . , Hence 
m£„(0, . . . ,0,t) vanishes on €o(C) for all ^ G Z. This implies that deg^ mi^ = D. As p is a monic 
polynomial of degree D in t, then p = m^^. So Iq is a primitive element of the extension K ^ L, 
and p is its minimal polynomial. 



Using the same notation of Lemma 3.1 we have 



^{U,L{U,x))x, = -^{U,L{U,x)) e C. 

As this identity only involves polynomials in Q[Uo, . . . , Ur][xi, . . . ,x„], it can be directly evaluated 
to obtain the parametrization w. □ 

In particular this shows that the total degree of the polynomials in the geometric resolution is 
bounded by degp < D and degWi < D (See also |27[ Prop. 3]). 

Lemma |3.2| can be directly applied to compute a geometric resolution of an equidimensional variety 



V which satisfies Assumption 1.2 from a given Chow form of V 



Corollary 3.3 Let notations and assumptions be as in Lemma 3.i. Suppose that there is given a 
Chow form Ty of V , encoded by a sip of length L. Then, there is an algorithm which computes a 
geometric resolution of V associated to t within complexity 0{nL). All polynomials arising in this 
geometric resolution are encoded by sip's of length 0{L). □ 

Lemma |3.2| also yields, from C/iy , a geometric resolution of the fiber Z associated to an affine linear 
form as ChziUo) = Chv{Uo, ei, . . . , e^). This is summarized in Procedure GeomRes (Subroutine 

In Procedure GeomRes (Subroutine |^), as we do in all zero-dimensional situations, we use the 
definition of geometric resolution stated in Section || to avoid divisions by p'. In line || of this 
subroutine, Res(/, g, ^1,^2) is a procedure that, using basic hnear algebra, computes (p, <?i, 92) where 
p is the resultant between the univariate polynomials / and g of degrees di and d2 respectively, 
and qi and 52 are polynomials of degrees bounded by 1^2 — 1 and di — 1 respectively satisfying 
p = qif + q2g. In line ^ Mod(/, g, c?i, (^2) is a procedure that computes the remainder of the 
division of the polynomial / of degree bounded by di by the polynomial g of degree bounded by d2 ■ 
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Subroutine 8: Computing a geometric resolution of a fiber 



procedure GeomRes{n,r, D,Chv,£,,c) 
^ n is the number of variables, 

^ r,D are the dimension and an upper bound for the degree of V respectively, 

# Chy is the normalized Chow form of V, 

# ^ := (^1, . . . , e A"^ such that = degV, where Z^:^Vn V{xi ~^i,...,Xr~ Cr)- 

# Co G s.t. £q :— Lq{cq,x) is the considered afhne linear form. 

# The procedure returns {Dq,p,v), where Do is the degree of V and {p,v) £ is the 

# geometric resolution of associated to £0 in case £0 separates the points in Z^. 
^ Otherwise, it returns error. 

1. P{UQ,t) := ChviiUoo - t, Uqi, . . . , J7o„),ei - ^leo, . . . ,6^ - ^rCo); 

2. (po,...,Pd) :=Expand(P(co,t),i,0,D); 

3. Do := D; 

4. while puf, = and Dq > do 

5. Do := Do - 1; 

6. od; 

7. {~l)^oP{co,t); 

8- (/o, 91,92) := Res(p,p',L'o,£'o - 1); 
9. if p = then 

10. return ("error"); 

11. else 

12. K, . . . , z«„) := ((-l)^°+iaP/9f/oi(co, t), . . . , (-l)^''+i9P/9f/o„(co, t)); 

13. ivi,...,vn) := (Mod(i 92^1, p,27:>o- 1,150),..., Mod(i q2Wn, p ,2Do - 1, Do)) ; 

return(i:»o,p, t^i, . . ■,Vn); 
end. 



Proposition 3.4 Let V C A" 6e an equidimensional variety of dimension r and degree bounded 
by D. Let (^1, . . . ,^r) G &e such that Z^ := n V{xi — . . . , Xr — ^r) is a 0-dimensional 
variety of cardinality deg V . Let be given Chy encoded by a sip of length L and the coefficients of an 
affine linear form £0 G Q[a;i, . . . ,Xn] which separates the points in Z^. Then, Procedure GeomRes 
( Subroutine ^ computes a geometric resolution of Z^ ( in the sense of Section |^ within complexity 
0{nD^L + D*). □ 



On the other hand, next result shows the converse of Corollary 3.3: To derive a Chow form from a 
given geometric resolution is quite standard in the zero-dimensional case, but it was by no means 
clear up to now how to generalize that for varieties of arbitrary dimension. Here we show how to do 
that within polynomial complexity. This is done by deriving, from the given geometric resolution of 
V, a geometric resolution of the fiber Z and a system of local equations for V at Z, which enables 
us to apply Procedure ChowForm (Subroutine 0). 

Proposition 3.5 Let V C A" be an equidimensional variety of dimension r and degree D which 



satisfies Assumption l.i. Let £0 £ Q[a;i, . . . ,Xn] be a linear form which separates the points of Z . 
Suppose that there is given a geometric resolution {p,w) of V associated to £ :~ {£o,Xi, . . . ,Xr), 
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encoded by sip's of length L. Then there is a bounded probability algorithm which computes (a sip 
for) Chv within complexity 0{n^^D^^{D + L)). 

Proof.- First we derive a geometric resolution of Z associated to ig: 
We know that Chz{Uo) = ChyiUo, ei, . . . , e^). Thus, 

rz{Uo)it) ^rv{Uo,ei, . . . ,er){t,0, . . . ,0) eQiUoM. 



The geometric resolution {p,w) of V associated to i is given by Lemma 3.2. Applying the same 
lemma to Z, we deduce that the geometric resolution (pq,wo) of Z associated to is po(t) := 
Pit, 0, . . . , 0) e Q[t] and wo{t) w{t, 0, . . . , 0) G Q[t]". 

Now, let us derive a system of local equations of V at Z: 

Let Ci £ Q"^^, r + 1 < i < n, be such that the afSne linear forms £i := + CnXi + • • • + Ci„x„ G 

linearly independent and such that each of them separates the points of Z. 
For r + 1 < i < n define 

H, := \p'{Mp)t- {p'i,{w/p')){Mp) 
where Mp e Q[xi, . . . ,Xr]^'^'^ denotes the companion matrix of p. Since p'ii{w/p') belongs to 

Q[xi,...,Xr][t], H,eQ[xi,...,Xr][t]. 

Observe that Xi — {wi/p'){£Q{x)) in L implies that in L 

= Qop'(^o) + C^lWl{£o) + ■■■+ C,nWn{£o) ^ (jl h(w j pJ . 

Thus, as My is the matrix of multiplication by with respect to K ^ i, we conclude that 

— Ip'(-^p) I where m^^ is the minimal polynomial of £i over K. 
The assumptions that li separates the points of Z and the projection tt : — > A'' , a: i-^ (x\, . . . , Xr) 
is finite, imply that m^. belongs to Q[a:i, . . . 
Therefore, for r + 1 < i < n, we can define 

:= rfmXk) = I ,/l.r ., Hi{xi, . . .,Xr)i£i) 
\p'{Mp)\ 

These are squarefree polynomials in separated variables which vanish over V, and so, it is easy to 
verify from the Jacobian criterion that fr+i, . . . , fn is a system of reduced local equations of V at 
Z. 

Observe that as there exist a, 6 G Q[a;i, . . . , Xr][t] such that discr(p) — a{t) p{t)+b{t) p' (t) , discr(p) Id = 
h{Mp)p'{Mp). On the other hand, as degZ = degF, discr(p)(0, . . . , 0) = discr(po) 7^ 0. Therefore, 
I h{Mp)p'{Mp) 1(0, ... ,0) = (discr(po))^ ^ 0. We conclude that \p'{Mp) |(0, . . . , 0) and hence, 
we can use the point (0, . . . , 0) to perform Procedure PolynomialDivision (Subroutine ^ in order to 
obtain division free sip's for fr+i, ■ ■ ■ , fn- 

Finally, we apply procedure ChowForm (Subroutine to Z and {fr+i, ■•■,/«}■ 

Let us decide now the random choices in order to insure that the algorithm has an error probability 

bounded by 1/4: 

We need Cr+i,...,c„ G Q"+^ satisfying the stated conditions of independence and separability. 
These conditions are satisfied provided that 



p(cr+i) • • • p{c„) 



i^ij ) r+l<i<n 
l<j<n— r 



^0, 



where p :— discrj Vz G Q[C/o] \{0}. As Vz is an homogeneous polynomial of degree D and degj Vz — 
D, degp < D {2D — 1). Thus the degree of the polynomial giving bad choices is bounded by 
{n — r) D {2D — 1) + {n — r). We choose i := 8n D^ in order to apply Schwartz lemma. 

Now we compute the complexity of the algorithm: 

The dense representation of the geometric resolution of Z associated to £0 is computed within 
complexity 0{nD^ L) (using Procedure Expand). 
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The construction of the random choice for the affine Hnear forms £r+i, . . . ,£n is not relevant here. 
The computation of each polynomial Hi requires 0{D^) operations for the computation of the 
determinant plus the computation of each coefficient of the matrix, that is 0{D^L) more operations, 
Hence, computing Hi requires 0{D^ {D + L)) operations. 

By Lemma 1.7, taking into account that the total degree of each fi is bounded by D (since it is the 
minimal polynomial of the affine linear form Hi), and that the lengths of Hi and |p'(Mp)| are of order 
0{D^{D + L)), the complexity of the final division for computing each f, is 0{D'^{D + D^{D + L))) = 
0{D^D + L)). 

Finally, Lemma U gives the final complexity ©(r* \og^{r D)n'^ D^'^ {D + L)) = 0{n^^D^^{D+L)). 
□ 



3.2 Intersection of a variety with a hypersurface 

Let V C A" be an equidimensional variety defined over Q, and let / S Q[a:;i, . . . , a;„] be a non-zero 
divisor modulo I{V). In this section, we compute, from the Chow form of V and the equation /, the 
Chow form of the set-theoretic intersection V fl V{f) C A". In order to do this, we use generalized 
Chow forms, which we define now. We refer to |43| and Section 2.1.1] for a more extensive 
treatement of these generalized Chow forms. 

We assume that dim = r and that deg/ < d. As before, for z = 0, . . . ,r, wc introduce a group 
Ui — (UiO, ■ • ■ J Uin) of n 4- 1 variables; we introduce also a group U{d)r of ('^'^") variables. We set 

Li := UiO + UtiXi H \- UinXn , -^r ^ U{d)ra 

\a\<d 

for the generic affine linear forms in n variables associated to Ui and the generic polynomial of degree 
d in n variables associated to U{d)r- 

Set N := r{n + l) + and \ei W d x V be the incidence variety of Xq, • ■ • , Lr-i,Fr with 

respect to V, that is 

W := {{uo, . . . ,Ur-i,u{d)r;0 e x A"; 

e e V, LoK, 0=0,..., L,„iK_i, = 0, Fr{u{d)r,0 = 0}. 

Let TT : A^ x A" denote the canonical projection onto the first coordinates. Then Tr{W) is a 

hypersurface in A^. A generalized Chow form or d-Chow form of V is any squarefree polynomial 
J'dy G Q[C/o, • • • , Ur-i, U(d)r] defining ^{W) C A^. 

A d-Chow form Tdy happens to be a multihomogeneous polynomial of degree d deg^ in each 
group of variables Ui, and of degree degT^ in the group U{d)r- 

If the variety V satisfies Assumption |1.2| , we define the normalized d-Chow form of V as the unique 
d-Chow form Ch^y € Q[Uo, . . . , Ur-i, U{d)r] of V satisfying Chdy{eo, ■ ■ ■ , Cr-i, eld)) = 1, where 
e{d) is the vector of coefficients of the polynomial xj?. 

Let V and V{f) denote the closure in P" of V and V{f) respectively. Set V n V{f) = Uc <^ 
for the irreducible decomposition of n V{f) C P" and, for each irreducible component C, let 
Tc G Q[?7o, ■ ■ • , Ur-^i] denote a Chow form of C. Then Q Prop. 2.4] states that 

Cha,v{Uo,...,Ur-i,f) = \XlT'^'' (7) 

c 

for some A £ Q* and some positive integers mc G N. (Here we wrote Chd.y{Uo, . . . , Ur-i, f) for the 
specialization of the group U{d)r into the coefficients of the polynomial /.) 
On the other hand, as F n V{f) — Uc!z:{a:o=o} ^^"^ polynomial 

n 

C^{xo=0} 
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is a Chow form of n V{f). 

Hence, in order to compute Chyf^yi^j-^, the goal is to compute first Chdy{f) :— Chdy{Uo, • ■ • , Ur-i, f) 
and then clean the multiplicities and the Chow forms of components contained in the hyperplane 
{xo = 0}. 

The following result enables us to compute a d-Chow form from the standard one. We recall some 
of the notation of Subsec tion 2.2: for an equidimensional variety V C A" of dimension r and degree 
D satisfying Assumption 1.2, we set K := Q{Uq, . . . , C/r-i) and I{V)° for the extension of the ideal 
of V to K[xi, . . . , Xn]- Also recall that 

V° V{I{Vr) n V{Lo, Lr^i) C A"(l?), 



is a 0-dimensional variety of degree D, and that Nyo refers to the Norm as defined in Section 2.1 

Lemma 3.6 Under Assumption |J.4 we have 

Chd.y ^Chv{Uo,...,Ur-i,eof f^vo{Fr). 

Proof.- Let Ch^yo G K[Ur] be the d-Chow form of V^. First one shows — exactly as in Lemma |2^ — 
that there exists G K* such that Chdy = XdCh^.v"- Set e{d)o for the vector of coefficients of 
the polynomial Xq. Evaluating this identity at U{d)o e{d)o we obtain 

ChdyiUo, . . . , Ur-i,e{d)o) = \dChdyo{e{d)n) = Xd- 

Consider the morphism Qd ■ Q[Uo, ■ ■ ■ , Ur-i, U{d)r] Q[Uo, . . . , Ur-i, Ur] defined by Qd{Li) = Li 
for < i < r — 1 and Qd{Fr) — Lf. Then gd{Chdy) — Chy (see Lemma 2.1]), which implies 
that 

Chdy{Ua, . . . , Ur^i,e{d)o) = Qd{Chdy){Uo, J/^-i, eg) = Chv{Uo, . . . , Ur-i,eo)'^. 

Therefore = Chv{Uo, . . . , C/, -i, 60)*^. The statement follows immediately from this identity and 
the observation that Chdyo = NYo{Fr). □ 



To clean the components of Vr\V{f) lying in the hyperplane {xq = 0} we use the following criterion: 

Lemma 3.7 Let W C P" be an irreducible variety of dimension r — 1. Then W C {xq = 0} if and 
only if Tw does not depend on the variable C/qo- 

Proof.- In case J-w does not depend on J/qo we have that 

Tw{eo, Ui,. . . , Ur-i) = 0, 

which is equivalent to the fact that W is contained in the hyperplane {xq = 0}. 

On the other hand, assume that W C {xq = 0} = P"~^. Then Tw coincides with the Chow form of 
W considered as a subvariety of this linear space, see e.g. the proof of Lemma 2.6]. Hence J-w 
does not depend on Uqo, and as a matter of fact, it does not depend on any of the variables Uio for 
0<i<r-l. □ 



Let again Tc G Q[Uo, ■ • • , C/r-i] denote a Chow form of an irreducible component C of VnV{f) C 
Recalling Identity 0, set 



^i-= n ^c"" and ^2:= n 

CC{xo=0} C<;t{xo=0} 



mc 

c ■ 



Then Chdv{f) ~ \T\Ti for A G Q* and the squarefree part (J-2)iod of Ti is a Chow form of 
V^V{f). ^ 

By the previous lemma, J-\ does not depend on [/goj while all the factors of Ti do. 
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Therefore 

dChayif) . ^ 
— A y"i - 



dUoo dUoo 

and so 



^ Chd,v{f) , , 



is a Chow form of n V{f). 



l.i. 



Lemma 3.8 Let V C A" be an equidimensional variety of degree D which satisfies Assumption 
and let f e Q[xi, . . . , a;„] of degree hounded by d he a non-zero divisor modulo I{V). Assume that 
Chy and f are encoded by sip's of length bounded by L. 

Then there is a bounded probability algorithm (Procedure Intersectfon (Subroutine^) below) which 
computes the Chow form Chy f^y(^j^ of the intersection variety VClV^f) within (worst-case) complexity 
0{{ndDf^L) . 



Proof.- Our first goal is to compute Chdy{f) S Q[Uo, . . . ,Ur-i] by means of Lemma |3.6| . To 
obtain Nyo(/) we derive first a geometric resolution of from its characteristic polynomial and 
Lemma It is easy to check that the polynomial 

pit) := i-lfChyiUo, . . . , Ur^i, [UrO " t, Url, . . . , Um)) 

is a characteristic polynomial of V'^ , with leading coefficient a :— Chy(Uo, . . . , Ur-i, eo). 
Then, the geometric resolution of associated to Lr is given by 

-p{t)eK[Ur][t] and -w{t) KprW where w := - ( ^ , . . . , ^ 

For 7 e F", if we denote by f^ the homogeneization up to degree d oi f with respect to a new 
variable a;o and p' the derivative of p with respect to t, we have 

p'{Lrh)r /(7) - f'\p'{LAj)), W,{Lrm, W^iLrh))). 

Thus, if M denotes the companion matrix of ^p(t), we get 

\p'{M)\'Nyoif) = \f'^p'{M),W,{M),...,Wn{M))\. 

In order to avoid divisions (since M e K[Ur]), we replace M by Mp := aM and p' , Wi, . . . , w„ by 
their homogeneizations (p')^, w^, . . . , up to degree D such that Mq :— a^p'{M) = (p')''(ald, Mp) 
and for 1 < i < n, Mi :— a^Wi{M) — wf(ald, Afp). Therefore, multiplying both sides by a''^ = 
ja''^ Idj'', we obtain 

|A/o|'Nyo(/) = |/"(Mo,Mi,...,M„)|. 



Finally, from Lemma 3.6, we conclude that 



ChMf) = a^^yo{f)=''^^l^^^^^^^ eQ[[/o,...,CZ.-i]. (9) 

We compute this quotient applying Procedure PolynomialDivision (Subroutine |^). 

Next we apply Identity |^ to compute a Chow form J- J-yni^(/) from Chd.y{f)'. we first compute 

the polynomial 

G gcdiChd,yif),dChd,y{f)/dUoo) (10) 

applying Procedure CCD (Subroutine ^ and then perform the division J- = Chd,y{f)/G applying 
again Procedure PolynomialDivision. 



Finally, as Assumption 1.2 holds, J-{eQ,...,er) ^ and we obtain the normalized Chow form 



Chy^V(^j:-) = T I T{eo, . . . , Cr) 
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Now let us check the size of the sets of points we have to take to insure that the algorithm has an 
error probability bounded by 1/4: 

First, in order to compute Chayif) we need u e Q('~+i)("+i) such that |A/o|(u) 7^ 0. But let 
us observe that in fact \Mo\{ei, . . . ,er,Ur) 7^ G Q[?7r] (so it is enough to choose randomly 
Ur & Q("+^) such that |Mo|(ei, . . . , e^, Mr) ^ 0). This is due to the fact that a{ei, . . . ,er) — 



Chv{ei, . . . ,er, cq) = ±1, thus Assumption 1.2 implies that C/iy (ei, . . . , e^, Ur) — ±Chz{Ur)- Hence, 
Pz{Ur, t) '■— p{ei, . . . , er){Ur, t) IS & characteristic polynomial of Z, whose discriminant docs not van- 
ish, and then, the polynomial |Mo|(ei, . . . , e^, Ur) ^ S Q[C/r]- 

Now, as deg |Mo|(ei, . . . ,er, Ur) < D^, if we take Ur := Random(n+l, 12 D^) we infer that with prob- 
ability at least 1 — 1/12, (ei, . . . , e,., u,.) is a good base point to apply Procedure PolynomialDivision 
and obtain Chdyif)- 



Next we compute G applying Procedure GCD [6 (1 -I- log 12)] = 26 times (see Remark 1.6 so that 
its error probabiliy is at most 1/12). 

Finally as G is a polynomial of degree bounded by rdD in r {n + 1) variables, choosing u :— 
Random(r (rt -|- 1), 12 rdD) we also guarantee that the probability that u is a good base point to 
perfom the last division is at least 1 — 1/12. 

Thus, the error probability of the whole algorithm is at most 1/4. 
Now let us compute the (worst-case) complexity of this algorithm: 

The whole complexity of computing the numerator and denominator in Identity is of order 



0{n{(fD^L + Z?")) = 0{n(fD'^L). By Lemma |L7| the complexity of computing Chd,v{f) is of 
order 0{{ndD)^{ndD + nd'^D^)L) = 0{n^d'^D'^ L). 

Then, we apply Lemma and Proposition |l.5| to compute a sip of length 0{'n7 d^D^'^L) for G of 
Identity within complexity 0{n^d^D''^L). 

Finally, when we perform the last division, the overall complexity of computing C/iyny(/) is of order 
0((ndD)i2L). □ 

We summarize the algorithm in Procedure Intersection (Subroutine 
3.3 Separation of varieties 

Let V C A" be an equidimensional variety of dimension r. Let g € Q[xi, . . . , Xn] \ {0}, and set Y for 
the union of the irreducible components of V contained in V{g) and W for the union of the other 
ones. Hence Y and W are equidimensional varieties of dimension r such that V — YUW, Y C V{g) 
and g is not a zero divisor modulo I{W). 

The following procedure (Subroutine |lO|) computes the Chow forms of Y and W from a Chow form 
of V and the polynomial g. 



For the sake of simplicity we assume that V — and therefore Y and W — satisfies Assumption |L 
Lemma 3.9 Let V C A" be an equidimensional variety of degree bounded by D which satisfies 



Assumption LI. Let g G Q[a;i, . . . , x„] \ {0} of degree bounded by d and Y and W defined as above. 
Assume that Chy and g are encoded by sip's of length bounded by L. 

Then there is a bounded probability algorithm (Procedure Sep ( Subroutine below) which computes 
the Chow forms Chy and Chw within (worst-case) complexity 0{{ndD)^L). 

Proof.- Let Vv S Q[C^o, ■ • ■ j . . . ,Tr] be the normalized characteristic polynomial of V , as 

defined in Subsection p.l| and set P' := dVy /dT^. 



We consider the following map, already introduced in Subsection 3.1 



(^y : a('^+i)("+i) X y ^ ViVv) , (^^0, ....ur;0^ ("o, ■ ■ ■ .u^, LqIuo, 0, ■ • • , Lr{ur, 0). 



By Lemma 3.1 tpv is a birrational map which in fact is an isomorphism when restricted to 
U := (A('^+i)(»+i) X V) \ V{P'{Lo, Lrj) ^ W :- V{Vv) \ V(P'), 
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Subroutine 9: Intersection with a hypersurface 



procedure Intersection(n, r, D, f, d, Chy) 

# n is the number of variables, 

^ r,D are the dimension and the degree of V respectively, 

# / G Q[a;i, . . . ,Xn] is a non-zero divisor modulo I{V) of degree bounded by d, 

# Chv is the normalized Chow form of V, 

# The procedure returns the normalized Chow form Chvnv{f) of the intersection variety VnV{f). 



4. Mp := a CompanionMatrix(p/a); 

5. for i from 1 to n do 

6. := Homog(wi,£>); 

7. Mi :=wf(a,Mp); 

8. od; 

9. ip'Y^ :=Homog(f 

10. Mo := (p')'^(a,Mp); 

11. /^:=Homog(/,d); 

12. Mf :=/''(Mo,Mi,...,M„); 

13. Hi :=\Mf\; 

14. H2 := |Mo|; 

15. Ur := Random(n + 1, 12 D^); 

16. if i/2(ei, . . . , e^., Ur) = then 

17. return( "error"); 

18. else 

19. Chd,v{f) := PolynomialDivision(a''iJi, rdD, (d, Cr, Mr)); 

20. G := GCD(CVv(/), dChd,vif)/dUoo, (Uo, . • . , C/r-i), rdl?; 12); 

21. u:=Random((r + l)(n+l),12rd£)); 

22. if G{u) = then 

23. return( "error"); 

24. else 

25. T := PolynomialDivision(G, Chii,v{f), rDd, u); 

26. Chvnv(f) ■= ^/^{eo, ■ ■■ , er); 

return(C/iyny(/)); 
end. 



with inverse 



2. a 



1. p 



{ — l)^Chv{Uo, . . . , Ur-l, {UrO — t, Url, . . . , Urn))', 

ChviUo,...,Ur-i,eo); 




t/jv ■■ {uo, ■ ■ ■ 



Ur; to,. . . tr) !-»• 



Uo, ...,Ur; - — 



p' da 



1 dVv 



P' dUon 



1 dVv 



) 
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Define 



G:= {Prrv{9)=9''{P': 



dVy 

dUol' 



dVv_\ 



where g'' :— Homog(5,d). Thus induces an isomorphism between V{g) D U and V{G) DU. 
Hence V{Vy) equals the union of the components in V{Vv) which are contained in the hypersurface 
V{G) C a('^+i)("+i)+('^+i), and ViVw) is the union of the other ones. As Vy is a squarefree 
polynomial we conclude that 



Vy :=gcd(G,n 



Pi 



w 



and therefore, from Identity (H) of Section 3.1, we obtain that 



Ty = Py(C/)(0) and Tw = Pw{U){0) 
are Chow forms of Y and W respectively. 

Note that as Vy \ Pv, 'Py (eo, . . . , e^, 0, . . . , 0) 0, thus e := (ep, . . . , e^, 0, . . . , 0) is a good base 
point to apply Procedure PolynomialDivision. Thus the only probability step of this algorithm is 
the computation of the Greatest Common Divisor between Vv and G. 

Now we estimate the (worst-case) complexity of the algorithm: 

The characteristic polynomial Vv can be computed from Chy with complexity 0{L) using Iden- 
tity (||) in Section Its partial derivatives with respect to Tq and C/qi, . . . , Uon can be computed 
within complexity 0{nL). 

The polynomial G is obtained within complexity 0{(P{d + nL)). 

As degVv = ('''+1) D and degG < d ((r+l)_D — 1), both bounded by (r+l) dD, the gcd computation 
of ■Py requires (ndD)^{d + L) additional arithmetic operations. 

^From Lemma 1.7, the polynomial division for Vw is then performed within complexity 0{{ndD)^ L). 



The final specialization T i-^ does not change this estimate. 

Therefore, the (worst-case) complexity of the algorithm is of order 0{{ndD)^L) 



□ 



We summarize the algorithm in Procedure Sep (Subroutine 10 



Remark 3.10 In case the variety V does not satisfy Assumption l.i, this procedure can be modified 
within the same hounds of complexity so that, from a Chow form of V , we obtain Chow forms of W 
and Y . The only problem that may appear in the previous lemma is that Vy{e) may be zero and we 
will not be able to accomplish the polynomial division. To solve this, we can modify Subroutine ^ in 
the following way: we choose a random point so that we can apply the polynomial division subroutine 
with error probability hounded by i and we change the error probability of the GCD computation also 
by I (by repeating it several times) in order that the error probability of the whole procedure is still 
bounded by j. 



3.4 Equations in general position 



The algorithm we construct in Subsection 3.5 works under some genericity hypotheses on the input 
polynomial system. This is one of the main reasons — but not the only one — for the introduction of 
non-determinism in our algorithm: there are no known efhcient deterministic procedures to obtain 
these hypotheses from a given system. In order to achieve them we replace the system and the 
variables by random linear combinations. Effective versions of Bertini's and Noether normalization 
theorems enable us to estimate the probability of success of this preprocessing. 

The complexity of our algorithm is controlled by the geometric degree of the input system, that is the 
maximum degree of the varieties successively cut out by the equations obtained by this preprocessing. 
To define this parameter, which is a suitable generalization of the geometric degree of a 0-dimensional 
system introduced in pj] , we first give the following definition: 
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Subroutine 10: Separation of varieties 



procedure Sep(n, r, D, g, d, Chy) 

# n is the number of variables, 

# r,D are the dimension and an upper bound for the degree of V respectively, 

# <?GQ[xi,...,a;„]\{0}, 

# is a bound for the degree of g, 

# Chv is the normalized Chow form of V. 

# The procedure returns the normalized Chow forms Chy and Chw- 

1. Vv '■= Chv{{Uoo — To, Uoi • • • , Uor), ■ ■ ■ , {UrO — Tr, Url, ■ ■ ■ , Urn))', 

2. 5^ :=Homog(5,(i); 

dVv 9Pv \ 



^' ^ ' ^ > ' ' QJJ^^ ' • • • ' QJJ^^ J ' 

4. Vy ■.= GCB{G,Tv,iUo,...,Ur,To,...,Tr),{r + l)dD); 

5. Vw '■= PolynomialDivision('Pi-, Vy, {r + 1) D, (eo, . . . , e^-; 0, . . . , 0)); 

6. Chy := Vy{U){0)/ry{eo, . . . , e^; . . . , 0); 

7. Chw := Vw{U){0)/Pw{eo, ■ ■ ■ ,er;0, . . . ,0); 
return(Cfty', Chw)', 

end. 



Definition 3.11 Let g G Q[xo, . . . ,ar„] be an homogeneous polynomial, Ig C Q[a;o . . . , a:„]g he an 
homogeneous ideal, and y C P" he a projective variety. We say that Ig is radical of dimension r 
outside Vg := V — V{g) if every primary component Q of Ig such that V{Q)g <^ Vg is prime of 
dimension r. 

An analogous definition holds for an ideal in Q[xi, . . . , Xn]g and an affine variety C A^. 

Let fi,...,fs,g € Q[xo,...,Xn\ be homogeneous polynomials of degree bounded by d, and set 
Vg = V(/i, . . . ,fs)g C ¥g. We assume that Vg ^ P^, that is fj ^ for some j. We also assume 
w.l.o.g degfj = d for every j: if this were not the case, we replace the input system by 

^d-deg/,-^^ , 0<i<n,l<j<s. 
For ai, . . . , a„+i e we set 

Qi{ai) := an fi-\ h Oisfs 

for the associated linear combination of /i , . . . , /s , which — by the assumption that deg fj = d — is 

also a system of homogeneous polynomials of degree d. 

Let A be the set of (n + 1) x s-matrices A = (oi, . . . ,a„+i 

Ii{A) := (Qi(ai), . . . ,Qi(aj)) C Q[xo, . . .,Xn], I <i <n+l, satisfy: 

. V{In+i{Aj)g = VginP^. 

• For 1 < i < n, if V{Ii{Aj)g ^ Vg, then Ii{A)g is a radical ideal of dimension n — i outside Vg. 

These are the first genericity hypotheses the polynomials should verify in order that our algorithm 
works. 

For every A G A we set S{A) := max{deg V{Ii{A))g ; 1 < i < n + 1}. 
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Definition 3.12 Keeping these notations, the geometric degree of the system 



A = 0, = 0, g^O 

is defined as 

S := Sih, ...,/.; 5) := max{S{A) ; A e A}. 
Note that Bezout inequality implies S < d". 

Remark 3.13 For a system of polynomials Fi, . . . , Fs,G €z Q[xi, . . . , Xn] (non-necessarily homoge- 
neous ) of degree bounded by d, the affine analogue (5aff of the geometric degree is defined in exactly 
the same manner, but without preparing the polynomials to make their degrees coincide. 
In fact, if for 1 < i < s, di := deg Fi, d := max^ di and Ff , G Q[a;o, • ■ ■ , a;„] are the homogeniza- 
tions of Fi and G respectively, then 

6,s{Fi, . . . , F,; G) = S{x^^-''^ F^ . . . , x^^^ F^ ; xo G'^). 

Let 

Vg = VaU---U Vn-i 

be the equidimensional decomposition of Vg in P^, where Vi is either empty or of dimension i, and 
let A = (ai, . . . ,a„+i)* G A. 

For i = 1, . . . as h{A) C (/i, . . . ,/^), Vg C V{h{A))g always holds. Moreover, i{V{h{A))g = 

Vg for some i, then V{Ij{A))g = Vg for all j > i. Also, observe that the ideal Ii{A) is generated by 
i polynomials, so every irreducible component of V{Ii{A)) has dimension at least n — i. Thus, we 
infer that for r := n — i, < r < n — 1, we have 

V{InMA))g = U l^, U • • • U 14-1 

where V^ is an equidimensional variety of dimension r. (We set V^ = for every r such that 
V{In~r{A))g ~ Vg since in these cases Vg = U • • • U Vn-) 
From now on, Qi(ai) will be denoted simply by Qi. 

The condition that A G A implies that, in case V^ ^ 0, Qn-r+i is not a zero divisor modulo I{V^)g. 
In this case, we have 

v;'_iuyr„iUv;u---uK-i - v{Qi,...,Qn-r+i)g 

= (K' u u • • • u v;_i) n v{Qn-r+i) 

- {Vlr\V{Qn-r+l))i^Vri^■■■i^Vn-l, 

as for all i, Vi C V{Qn-r+i)- Hence, since dim(V^' n V{Qn-r+i)) = r — 1, we deduce that 

K n V{Qn-r+l) = K'-l U Vr-l U Vr-l (11) 

where 

V^-i = |J{ G; G component of n F((3„_r+i) n (14 U • • • U K-i) of dimension r - 1 } 
is an equidimensional subvariety of 14 U- • -UKi-i of dimension r — 1. We set 14- 1 := and I4i 
Now for bo, . . . ,bn G Q"+^ we consider the linear change of variables 

Vkibk) ■=bkoXo-\ i-bknXn , < fc < n. 

We say that {bo, . . . , 6„) is admissible if, under this linear change of variables, for < r < n — 1, 



• the varieties I4' U 14 U Vr satisfy Assumption 1.2 
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• the polynomials Qi, . . . , Qn~r € ^(KOg ^^'^ ^ system of local equations of at Zr := fl 
V{yi, . . .,yr). 



We construct the polynomials Qi, . . . , Qn+i and the variables yo, ...,?/„ by choosing the coefficient 
vectors a^, 1 < i < n + 1, and 6^, < fc < n, at random in a given set. In what follows we estimate 
the error probability of this procedure: 

Lemma 3.14 Let notation be as in the previous paragraphs. Let N he a positive integer and let 
fli g [0,87V(d+ 1)2")'' l<i<n+l, , 6fc e [0,27Vn2rf2»)"+i < /c < n 

be chosen at random. Then the error probability of A := (oi, . . . , an+iY being in A and (bo, . . . , 6„) 
being admissible is bounded by 1/N. 

Proof.- The set of matrices A contains in fact a non-empty open set of q("+i)^*: by the effective 
Bertini's theorem in Lemmas 1 and 2] or a local version of ||3^, Prop. 4.3 and Cor. 4.4] there 
is a non zero polynomial F with degi^ < 4 (d + 1)^" such that F(ai, . . . ,an+i) ^ implies that 
A = (ai, . . . ,a„+i)* e A. 

Assume now that ^ € A. By the effective Noether theorem version of |38, Prop. 4.5] there is a non 
zero polynomial G G A:[C/o, . . . , C/„] with 

n-l 

deg G < 2 ^ r deg(T/; U K U Vrf 

T'=0 

such that G{bo, . . . , 6„) =/= implies that under the linear change of variables given by (69, . . . , 6„), 
the varieties Vr U U Vr satisfy Assumption 1.2. Since, from Identitiy [ll|, 

deg(V;,' UVrUVr)<d degV;+i < d", 

degG < n(n- 1)^^". 

Now we will define a polynomial H € A:[?7i, . . . , /7„--i] such that H{bi, . . . , 6„_i) 7^ implies that 
the second condition for admissibility holds. 

Fix r, < r < n — 1. We know that {Qi, . ■ . , Qn-r)g is a radical ideal of dimension r outside Vg 
whose associated variety coincides with ]/' outside Vg. Thus, locahzing at any ^ G V^, ^ ^ Vg, we 
get ((Qi, . . . , Qn-r)g){ = J ^hat is, Qi, . . . , Qn-r is a system of local equations of V^ at ^. 

Therefore, it suffices to take new variables yo, ■ ■ ■ ,yn such that VJ^ CiVgCi V{yi, . . . , y^) = 0- 
From the definition of Vq, it is clear that Vq DVg = 

For l<r<ri — 1, asV^is definable by polynomials of degrees bounded by d and no irreducible 
component of V^ is contained in Vg, there exists a polynomial gr G k[xi, . . . ,Xn\ with deg((/r) < 
such that Vg C V^(.gr) and V^. fl V^(.g-r) is equidimensional of dimension r — 1. Let G fc[[/i, . . . , C/^] 
be a Chow form of VJ,' fl V{gr). 

Set iJ := Ylr=i -^r ^ ^[C^ii ■ • • J C^n-i]- The condition -ff (61, . . . , 6n-i) 7^ implies that, for every 
1 < r < n — 1, Tr{b\, . . . , 6^) ^ and so, 

Vl ^Vg(^ F(yi(6i), . . . , yr(br)) C K' n V{g,) D V{yi{bi), . . . ,yr{br)) = 0. 
Observe that H is a, non zero polynomial with 

n—1 n—1 n—1 n — 1 

deg = ^ deg :Fr^Yl n V{gr) < ^ d deg < XI ^""'"^^ < ('^ - l)^''- 

r— 1 r— 1 r— 1 r— 1 

Therefore, there exists a non zero polynomial condition of degree bounded by n{n—l)d'^"+{n—l)d^ < 
n'^d'^^ which ensures that the matrix {bo, ... , 6„) is admissible. 

The conclusion follows as usual from the Zippel-Schwartz test. □ 
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3.5 Proof of Theorem 1 



Let fi,...,fs,g e Q[xQ,...,Xn\ be homogeneous polynomials of degree bounded by d, and set 
Vg :— V{fi, . . . , fs)g C Set S for the geometric degree of the system fi = 0, . . . , fs ~ 0, g ^ 0. 

The algorithm is iterative and consists of two main steps, besides the preparation of the input 
(equations and variables). 

The preparation of the input enables us to work with an affine variety W instead of the input quasi- 
projective variety Vg and local systems of equations of certain auxiliary varieties appearing in the 
process. 

The first main step computes recursively the Chow forms of a non-minimal equidimensional decom- 
position of W. Here the crucial point which controls the explosion of the complexity is that the size 
of the input of an iteration does not depend on the size of the output of the previous step: the input 
of each recursive step has the same controlled size. 

The second main step clears extra components and computes the Chow forms of the equidimensional 
components of the minimal decomposition of W from which the Chow forms of the equidimensional 
components of Vg are obtained straightforwardly. 

This is a bounded error probability algorithm whose expected complexity is of order s{nd5)'~^^^^ L. 
Its worst-case complexity is s{nd'''^)^'^^^ L. 

For the rest of this proof, we set N := d^^". 

Input Preparation 

Set Vg = Vq\J ■ ■ ■ yjVn for the minimal equidimensional decompositon of Vg , where each Vr is either 
empty or of pure dimension r. 



First, applying Procedure Deg described at the end of Section 1.3 to fi, ■ ■ ■ , fs, we compute with 
error probability bounded by 1 /67V the exact degree of the polynomials fi, ■ ■ ■ , fs within complexity 
0{sd?L + n\og{sdN)). This also states whether these polynomials are the zero polynomial and, 
therefore, whether Vg = P^. In that case J^y^ = |([/o, ■ • ■ , Un)\ and for i < n, — 1. 
Thus, with error probability bounded by 1/QN we can assume we know the exact degree of the 
polynomials /i, . . . , /s, and that T/„ = and dim Vg < n — 1. 

We consider the polynomials 

hj--=xi~''"^^' f, , 0<t<n,l<j<s, 
hence we have now t < {n + l)s polynomials fij of degree d, that we rename /i, . . . , /j. 



We apply Lemma 3.14 to choose randomly a matrix A = (oi, . . . , a„+i)* G and a matrix 

B = {bo, . . . , e qI^+i)x(«+i) such that the error probability that ^ G A and B is admissible is 
bounded by 1/6A^. 

We can assume thus that the linear combinations (Qi, . . . , Qn+i) = A{fi, . . . , ft) satisfy 



Vg=V{Ql,...,Qn+l)g 

and, for < r < n — 1, 

• (Qi, . . . , Qn-r)g IS either empty outside Vg or a radical ideal of dimension r outside Vg. 

• y{Qi, ■ ■ ■ 1 Qn-r)g = VJ,' U U • ■ • U 14,-1, whcrc V^ is either empty or an equidimensional 
variety of dimension r with no irreducible component included in U • • • U Vn~i- 

• V^ n V{Qn-r+i)g = K-'-i ^ U Vr-1, where Vr-1 is either empty or an equidimensional 
variety of dimension r — 1 included in 14 U • • • U We set V^ := to extend this property 
to r = n. 



We can assume moreover that the change of coordinates y = B x verifies 
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• BV^ U BVr U BVr satisfies Assumption ^ 

• Qi{B~^y), . . . , Qn-r{B~^y) is a system of local equations of BVJ^ at BV^ fl ^^(yi, . . . , j/r)- 

The complexity of constructing the random matrices A and B and the inverse of the matrix B is of 
order C(sn''(log7V + logd)). 

Now, Assumption 1.2 implies that the varieties have no irreducible component at infinity. Hence we 
restrict to the afHne space: we set yo — 1 and denote by gi, . . . , Qn+i, h the set of polynomials in the 
new variables obtained from Qi, . . . , Qn+i, 9, that is: 

(gi, ... ,(7„+i) = A (1,2/1, ... ,?/„)) , ft. = g(B"^(l, 2/1, ... ,?/„)), 

where F := (/i,...,/t). 
We define 



W := V{qi, . . . , qn+i)h = BVnAl(lA". 

Let W = Wq U ■ • ■ U Wn-i be the minimal equidimensional decomposition of W, where for < r < 
n — 1, Wr is either empty or of dimension r, and let and Wr defined by the same construction 
as and Vr before, that is 

• Viqi, qr,-r)h = U W^, U • ■ • U W^„_i 

. n V{q^-r+l)h = W^;_i U W^,_i U Wr^i 

As the identity 



Wr ^ BVr^Kl 



holds, from a Chow form of Wr we obtain a Chow form of the corresponding Vr by means of the 
change of variables: 

J-y,(C/o, ...,[/,)= ChwAUoB-\ . . . , C/,B-i). 



We observe that W^ = BV^ D AJJ, and then qi, . . . , qn-r is a system of local equations of W^ at 
M^;nF(2/i,...,2/.). 



The error probability of this preparation step is bounded by 1/3N. Once the matrices A and B are 
fixed, we have that the complexity of computing the polynomials gi, . . . , qn+i, h and their length are 
all of order 0{sn'^dL). 

First Main Step 

From r = n — 1 to 0, the algorithm computes the Chow form of Wr U Wr and a geometric resolution 
of the fiber Zr := W!^ fl V{yi, . . . ,yr) (which also gives the degree Dr of W^). The former will be 
the input of the second main step while the latter is the input of the next step in this recursion. 
Each step of this recursion is a bounded probability algorithm whose error probability is bounded 
by l/ZnN provided that the input of the iteration step was correct. 

We begin with the fiber Z„ ~ V{yi, . . . , ?/„) — (0, . . . , 0) and its geometric resolution {t, {t, . . . , t)) 
associated to ^ = a;i. We also set Z?„ 1. 

Now, we are going to describe a step of the recursion. From a geometric resolution of Zr+i we 
compute a Chow form for Wr U Wr and a geometric resolution of Zr, which is the input of the next 
recursive step. Set Dr+i for the given estimate of deg W/_,_i. 



• Computation of C/ivF'^j^ : 

From the geometric resolution {pr+i, (ui, . . . , associated to the affine linear form ir+i of 
Zr-i-i, and the system of local equations gi, . . . , qn-r~i of W^.-^ at Zr-t-i, we compute the Chow 
form of W^_|_i applying Procedure ChowForm (Subroutine H). This step of the algorithm is 
deterministic and computes Chw'^^ provided that the polynomials and variables satisfy the 
genericity conditions and that the geometric resolution of Zr+i is accurate. Observe that by 
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Main Lemma 2.3 applied to the local system of equations qi, . . . , qn-r-i of^ degree d and length 
0{s'n?dL), the complexity and the length of the output are both of order 

HChw;^^) - 0{{r + lflog^{{r + l)Dr+iWd^Dll,isn^dL)) ^ 0{sn\ndDr+if'' L) . 
Computation of C/ivi/'^^nv((j„_r): 

Now we apply sufficient times Procedure Intersection (Subroutine ^ to compute the Chow 
form of W^_|_i n V{qn-r) with error probability bounded by l/18niV: by Lemma 3.S, the 
length of the output Chow form and the complexity of one iteration are both of order 

L(C/i„.;^^ny(9„_„)) = 0{{ndDr+iY^L{Chw:.^J) = 0{sn%ndDr+if^L), 



while, from Corollary 1.6 for the choice s = [6(log(18nA^) + 1)], the complexity of this step is 
of order 

0(((r + l)(n + 1) + 1) log(18n7V)(L(C/iiv;^ J + L(C V,^^ny(g„_.))) + log'(18niV)) = 

= 0{\og\N)sfv>{ndDr+if'^L). 

Computation of Ch^^^^^^^^: 

Observe that each irreducible component of W^+i ^^{in-r) is either an irreducible component 
of W!^ U Wr U Wr or an irreducible variety included in V{h). Therefore, we apply sufficient 
times Procedure Sep (Subroutine 10) to compute the Chow form of W!^ U Wr U Wr with error 
probability bounded by l/18niV: by Lemma 3.9| , the length of the output Chow form and the 
complexity of one iteration are both of order 

^(C/iiy;uiy.uv^.) = 0{{nd{dDr+i)fL{Chw,^^^vi,^^.))) = 0{sn^d\ndDr+if^L), 
while the complexity of this step is of order 

O ( log^(7V) sn^ d^ {ndDr+if^ L) . 

Computation of Chwi^ and Ch^ : 

Next, since Wr U Wr C V{qn-r+i) and no component of W^ does, we use g„-r+i to separate 
Chyy^ from Ch^^, . We apply sufficient times Procedure Sep (Subroutine 10) to compute 

the Chow forms of W^. and Wr U Wr with error probability bounded by l/18nA^: the length 
of the output Chow forms and the complexity of one iteration are both of order 

L(C/jv^;,CV,uw.) = (^(M(^^-+i))'^(CV;uw.uH>.)) = Oisn^d'^indDr+i^L), 
while the complexity of this step is of order 

0{log\N) sn^d^^{ndDr+i)^'^L) . 

Computation of a geometric resolution of :— W^. r\V{yi, . . . ,yr): 

We apply here Procedure GeomRes (Subroutine |^). It requires a random choice of the coeffi- 
cients of a separating linear form ir- We do that in order to insure that the error probability 
is l/6nN. The condition that a linear form separates the points of the fiber Zr is given by a 
polynomial of degree bounded by (^''l^'') < — - as AegZr < d"^"^. So we choose the set 
of coefficients of ir in [0, SnA^d^'^"^''')"^^. The complexity of constructing these coefficients is 
thus of order 0{{n + l)(log(niV) + (n — r) logd)) = C'(n^(log + logd)) and the complexity 
of computing afterwards the geometric resolution of Zr (that is, all its constant coefficients) 
adds, as Dr < dDr+i, 

0{n{dDr+ifL{Chw;) + d^D^^i) = 0{srt' d^^ {ndDr+if^ L) 

operations. 
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Summarizing, from the geometric resolution of Zr+i and the polynomials qi, . . . , Qn-r, the algorithm 
produces, within complexity 0(^\og\N)sn'^d^^{ndDr+i)'^^L) , all the coefficients of the geometric 

resolution of Zr and a sip of length O (^sn^ d^^ {ndDr+iY'^ for the Chow form of Wr U Wr- The 
error probability that the computed objects are not the correct ones, provided that the input was 
right, is bounded by l/3nN. 

Therefore, provided that the input preparation was correct, this algorithm is expected to compute 
Ch^ UM> , for < r < n — 1, with error probability bounded by 1/3-/V, within complexity of order 

n-l 

0{\og\N)sn'd''^i J2 {ndDkf^)L), 

fe=r+l 

and, by the iterative character of the algorithm, to compute all Ch^ , < r < n — 1, within the 
same complexity than that of computing Ch^^^:^^. 

Second Main Step 

For 0<r<ri — 1, in order to extract from the Chow form Ch^ the factor C/iw^, we define 
a hypersurface V{Gr) such that, probabilistically, Wr is exactly the union of all the irreducible 
components of Wr U Wr contained in V{Gr), and then we apply Procedure Sep (Subroutine |l^) to 
compute C/ivKr- 

Fix fc, l<fc<n — 1. We define a polynomial H^. G . . . , ?/„] such that, with error probability 

bounded by l/6(n — the following conditions hold: 

1. W^fcUW^fe CT/(iJfe), 

2. no irreducible component of Wr is contained in V{Hk) for r = 0, . . . , fc — 1. 

Let V be the characteristic polynomial of Wk U Wt- For any affine linear form = -^o(co,2:), we 
have that := P(co, ei, . . . , efc)(^0: yii • ■ • j 2/fc) vanishes on W^ U Wk- We determine now randomly 
such that Condition 2 holds with error probability bounded by l/6(n — 1)A^. This is a standard 
argument that can be found for instance in |2^, Section 2.3.5]): 

For any irreducible component C of VFoU- • -UW^fc^i there exists := £,n) G C —{Wk^Wk)- 

Now, if a linear form satisfies that for any f £ {Wk^Wk)C\V [yi—S^i , . . . ,yk—^k) (which is a zero- 
dimensional variety of degree bounded hydS), 4(C) ^ io{£,c), 'P{co, ei, . . . , ek){io{^c),ii , ■■■,ik)¥' 
holds. Hence C is not included in V{Hk). 
The condition to be satisfied is thus given by 

where C runs over the irreducible components of Wo U • • • U Wk-i and ^ G {Wk U Wk) fl V{yi — 
Cf , . . . , j/fc — C^). The polynomial has degree bounded hy d6^ < since deg Wq U • • • U Wk-i < S. 

Choosing cq '■— (0, cqi, . . . , con) G [0,6{n — l)iVd^"+^)", the probability that Hk does not satisfy 
Condition 2 is bounded by l/6{n — 1)N. Therefore the probability that, for 1 < fc < n — 1, at least 
one Hk does not satisfy Condition 2 is bounded by 1/6N. 

Now, for r — 0, . . . , n — 2 we define Gr := Ilfe^r+i ^k- Clearly, as Wr C Wr+i U • • • U Wn-i, Gr 
vanishes on Wr by Condition 1. On the other hand, as, by Condition 2, no irreducible component 
of Wr is contained in V{Hk) for r + 1 < fc < n — 1, splits Wr and Wr- 

For 1 < fc < n — 1, deg Hk < dDk+i and L{Hk) = 0(L{Ch^^^^^)) since we derive V from 
the corresponding Chow form by Identity Hence L{Hk) ~ 0[sn^d^^{ndDk+i)'^'^ L) . Thus, for 
< r < n — 2, degGr < dJ2k>r+i ^k+i and 

L{Gr) =0{sn^d^'^{ ^ {ndDk+if")L). 

k>r+l 
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The computation of all ff^, 1 < fc < n— 1, involves the computation of the random coefficients of each 
linear form (q, that is 0{n'^ {log N + n\ogd)) operations for each one of them, plus the complexity of 
computing and specializing each characteristic polynomial. Thus the total complexity of computing 
all Hk is of order ©(n^logiV + sn^d^^(^j,>2("c?^fc)*")-^) • We conclude that the complexity of 
computing all Gr, < r < n — 2, is also of the same order. 

This algorithm is expected to compute the right polynomials Go, . . . , Gn~2, provided that the Input 
Preparation and the First Main Step were correct, with error probability bounded by 1/6N. 

Now we apply sufficient times Procedure Sep (Subroutine ^o|) to C/i^ and G,. in order to 
compute Ch\Y^ with error probability bounded by 1/QnN: the length of the output Chow forms and 
the complexity of one iteration are both of order 

n-l 

L{Chw^)^0{{n{d Dk+i){dDr+i)fL{Ch^^^^^,Gr))^0{sn'd^\ndDf^L), 

k=r+l 

where D = max{Dfe : 1 < A: < n — 1}, while the total complexity of computing all Chw^ with error 
probability bounded by 1/6N , provided that the polynomials Gq, . . . , G„_2 were correct, is of order 

0{ \og\N)sn^^d^^{ndDf^L + slog^{s) log(d)L). 

Thus, the total error probability of the second main step is bounded by 1/3 A''. 

Finally, the Chow form J-y^ is obtained by changing variables back. This computation does not 
change the order of complexity involved. 

The total error probability of the whole algorithm is bounded by Moreover, in case each 

of the random choices was right, < S for every k, and therefore the Chow forms Tvr- of the 
equidimensional components Vr of Vg are encoded by sip's of length 

L{TvJ = 0{srJd^^{nd5f^L), 

and computed within complexity 

0{ \og\N)sn^^d^^{nd5f^L + ) . 

Since in any case, Dk < d^^^ < d"^^ for every 1 < fc < n — 1, the worst-case complexity of the 
computation is of order 

0{\og\N)sn^''d^^d^^"'L). 
Therefore the expected complexity of the algorithm is 

0[{l - l^)(log\N)sn^^d^^{ndS)^^L) + ^{\og\N)3n^'^d^'^d^^''L)). 

Fixing N :— d^^" , we conclude that the expected complexity of our bounded probability algorithm 
is of order 

0{log^{d^^'')sn^^d^^{ndSf^L + log^ {d^^'')sn^'^d^^L) = s{ndSf^^^L, 
while the error probability is bounded by 

We summarize in Procedure Equidim (Subroutine |ll|) the algorithm underlying the Proof of Theorem 
1. 
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Subroutine 11: Equidimensional decomposition 



procedure Equidim(n, d, /i, • • • , /s, g, x) 

# fi, ■ ■ ■ , fs,g are homogeneous polynomials in Q[a;o, . . . , Xn] and x := {xq, ■ ■ ■ , a;„); 

# d is an upper bound for the degrees of /i, . . . , fs, g. 

# The procedure returns the Chow forms J^Vo ) • • • > -^y„ of the equidimensional components 

1. N := d^S"; 

2. {di,...,ds) := Dcgifi, . . . , fs, X, d;6sN); 

3. if (di, . . . , ds) := (—1, . . . , —1) then 

4. {J^vo,. ■ • := (1, . . . , 1, {{Uo, Un)\); 

5. else 

6- (-^0 ^/li-'-J-^n fli ■ ■ ■ ■> ^0 " fs-i ■ ■ ■ iXf^ fs)] 

7. A := RandomMatrix(n + 1, s(n + 1), 48iV(d + 1)^"); 

8. B := RandomMatrix(n + 1, n + 1, UNn'^d'^'"); 

9. {yo,...,yn):=B{xo,...,Xn); 

10. iqu...,qn+i) := AF(B-i(l,yi,...,y„)); 

11. h:=g{B-\l,yu...,yn)); 

12. := 1; 

13. (c("),i?„,p„,,/")) ■.= {e,,l,t,{t,...,t)); 

14. for i from 1 to n do 

15. r := n — i; 

16. Chw^_^_^ ■■= ChowForm(n, r + 1, Dr+i, c^'^~^'^\Pr+i, v^'^+'^'> ,qi, Qn-r-u d); 

17. := Intersection(n, r + 1, -Dr+i, 9n-r, Chw:^^^ ; ISnA''); 

18. ^V^uw.uw', •= (Sep(n,r + l,rf£)r+i,/i,rf, J^;lB,nN))2\ 

19. (C/iw-;,C/i^y.^LJW^,) := Sep(n,r,d£)r+i,d, Ch^^^^^^^j^^, qn-r+vA^nN); 

20. c('') := Random(n + 1, 3nA/'d2("-'-)); 

21. (D^,^^,?;^) := GeomRes(n,r,d£>r+i,C/iw-;,(0, ...,0),c('')); 

22. od; 

23. for k from to n — 1 do 

24. Vk ■= ^hy^^^^j^^{{Uoo - To, Uoi, Uon), {UkO - Tr, Ukl, Ukn)); 

25. uf'^) := Random(n + 1, 6(n - l)iVrf2"+i; 

26. Hk ■.= Vk{u'^'^\ei,...,ek){uf^ +uf\i H V u:h\n,yi, ■ ■ ■ ,yk); 

27. od; 

28. for r from to n — 2 do 

29. Gr := riLr+i 

30. C/iiy^ := Sep(n, r,dDr+i,Gr, d{Dr+2 -\ h -Dn), C/i^^|_|^^; 6niV); 

31. J^Vr-=ChwAUoB-\...,UrB-^y, 

32. od; 

33. fi; 
return(jrv-g,...,jry^); 

end. 
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4 Applications 



We present some algorithmical applications of our results, concerning the computation of resultants 
and the resolution of generic over-determined systems. 

4.1 Computation of resultants 
4.1.1 The classical d- resultant 

As a first application of our results, we compute a sip for the classical resultant of rt + 1 generic 
homogeneous polynomials of degree c? in 1 variables. The algorithm follows directly from Lemma 



2.3 and is therefore deterministic. For the definition and basic properties of the classical resultant 



we refer for instance to [|13|, Chapter 3]). 

Corollary 4.1 There is a deterministic algorithm which computes (a sip for) the classical resultant 
Res„ d of n + 1 generic homogeneous polynomials of degree d in n + 1 variables within complexity 

Proof.- It is a well-known fact that the resultant Res„_d is the Chow form of the Veronese variety 
V(n, d) defined as the image of the morphism 

where TV := - 1. We recall that K, d is an irreducible variety of dimension n and degree d". 

We compute here the resultant by defining a system of local equations at an adequate fiber of V{n, d) 
in order to apply Lemma |2.3[ 

Let {ija ■ a e Ng^^, |q;| = d} be a set of homogeneous coordinates of and consider the projection 

TT : V{n, d) P" , {ya)a ^ {Vdeo ■ ■ ■ ■ ■ ydej 

where Cj is as usual the {i + l)-vector of the canonical basis of Q""'"^. This projection is finite 
^ Ch. 1, Thm 5.3.7]. Moreover, Z 7r-i((l : 1 : • • • : 1)) verifies that Z = (pnA^o) with 
Zo := {{1 : LUi : ■ ■ ■ : w„); = 1 for 1 < i < n }. Thus #Z ^ d" ^ degV{n,'d) and the 
n-dimensional variety V{n,d) satisfies Assumption 1.2 for the fiber Z. 

Let us define now a system of local equations of Vn,d at Z: For every a — {ao, . . . , a„) G (No)"+-^ 
such that |a| = d and a ^ {d — l)eo + e.^ (0 < i < n) we consider the polynomial 

These are N ~ n non-zero homogeneous polynomials of degree d — ao which vanish at Vn,d since 

From the Jacobian criterion one also checks that, as — y^e^^ "° and = for /? 7^ a and 
/3 ^ (d — l)eo 4- Bi, the Jacobian matrix of the system has maximal rank A'^ — n at any ^ G Z. 
Observe that the equations fa can be encoded by sip's of length 0{d). 

Next step in order to apply Lemma ^.3| is to compute a geometric resolution of the fiber Z. For that 
aim we compute its characteristic polynomial (considering it as an affine variety in {ydeo 7^ 0}) and 



apply Lemma 3.2 for a separating linear form. 

Let L :— UaVa be a generic linear form in A' -I- 1 variables, and let P = X]|a|=ci ^q^;" be the 

generic homogeneous polynomial of degree d in n -I- 1 variables associated to L. 
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The characteristic polynomial of Z is 



Vz{u,T) = [](r-L([/,0) = n (r-i(f/,(/p„,d(i,w))) 

= n (7^-m(i,^))) = n 

where P°(;7,w) = P(C/, (l,w)). Therefore, if we set A := Q[xi, . . . ,a;„]/(a;f - 1, . . . ,a;f, - 1), is 
then computed as the characteristic polynomial of the linear map A ^ A defined hy g P'^g within 
complexity dP^'^\ 

Finally, an easy computation shows that the linear form £ = ydeo +'^y(d-i)eo+ei +' ' ■ + '^"2/(d-iiefl+e„ 
separates the points in Z. Thus £ yields a geometric resolution of Z and we apply Lemma 2.3 to 
compute Res„ d within the stated complexity. □ 



4.1.2 Sparse resultants 



Let A — {ao, ■ ■ ■ ,aN} C Z" be a finite set of integer vectors. We assume that Z" is generated by 
the differences of elements in A. 

For < i < n let [/i be a group of + 1 variables indexed by the elements of A, and set 

F,:=Y,U^ax'' eQ[C/,][xf\...,xf ] 

aeA 

for the generic Laurent polynomial with support in A. Let Wa C (P^)"+i x (C*)" be the incidence 
variety of Fq, . . . , F„ in (C*)", that is 

Wa = {{^yo, . . . , ^^n; e (P^)"+i x (C*)" : F,(i.„ = 0, < i < n}, 

and let tt : (P^)"+i x (C*)" (pW)"+i be the canonical projection. Then 7t{Wa) is an irreducible 
variety of codimension 1. The A-resultant Res^ is defined as the unique — up to a sign — irreducible 
polynomial in Z[J7o, . . . , Un] which defines this hypersurface (see ||l9|, Ch. 8, Prop.-Defn. 1.1]). 
This is a multihomogeneous polynomial of degree Vol(^) in each group of variables Ui, where Vol(^) 
denotes the (normalized) volume of the convex hull Conv(^), which is defined as n\ times its volume 
with respect to the Euclidean volume form of R". 

Consider the map 

(C*)"-^P^ , (^"° : ••• :^""). 

The Zariski closure of the image of this map is called the afftne toric variety Xa C P^ associated to 
A. This is an irreducible variety of dimension n and degree Vol(^). Its Chow form coincides — up 
to a scalar factor — with the sparse resultant Res^ G Z[[/o, ■ • ■ , Un] (see Ch. 8, Prop. 2.1] and 
H Ch. 7, Thm. 3.4]). 

For a broader background on toric varieties and sparse resultants we refer to and . 

We apply the algorithm underlying Theorem 1 to compute the sparse resultant Res^ for the case 
that A C (No)" and the elements 0, ei, . . . , e„ — that is the vertices of the standard simplex of IR" — 
lie in A. To do so, we construct a set of equations which define Xa in the open chart {P^)yg, where 
ivo : • • • : 2/Jv) is a system of homogeneous coordinates of P^, and compute a Chow form of this 
variety. 

Corollary 4.2 Let A C Nq be a finite set which contains {0, ei, . . . , e„}. Then there is a bounded 
probability algorithm which computes (a sip for) a scalar multiple of the A-resultant Res^ within 
(expected) complexity (n + Vol(^))'-'("'^^ . 
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Proof.- W.l.o.g. we assume that in A, ao ^ and ai = Ci for z = 1, . . . , n. Set d niaxQg_4|Q;|. 
For 71 + 1 < j < iV we set 

f, := y',-'y, - y^^-^^y'^^' ■ ■ ■ y^- G Q[yo, • . • , y^]. 

Then, XA\{yo = 0} = V := V(/„+i, . . . , /7v)yo C (P^)yo- Therefore the Chow form of coincides 
with the one of V and can be computed by apphcation of Procedure Equidim (Subroutine 11) to 
the polynomial system /„+i, . . . , /at; yo- 

Each polynomial fj,n + l<j<N, can be encoded by a sip of length 0{d). Moreover, as for 
each a £ A, \a\ = Vol({0, ei, . . . , e„, a}) < Vol(^) since {0, ei, . . . , e„, a} C ^, then d < Vol(^). 
Therefore L{fj) < 0(Vol(^)) for n + 1 < j < iV. 

Now, as the toric variety is non-degenerated (that is, it is not contained in any hyperplanc in 
P^), H, Cor. 18.12] implies that 

iV + 1 < dimXA + degX^ = n + Vol(^). 

This gives an estimation for the parameter N. 

Finally, we have to estimate the geometric degree (5(/„+i, . . . , /jv; J/o)- As we want to compute 
this degree outside {yo = 0} it is enough to deal with linear combinations of the dehomogeneized 
polynomials fj obtained by specializing yo = 1 in the original fj for n + 1 < j < N . 
For 1 <i < N, n + 1 < j < N and Uij £ Q we set 

Qi := ttin+i fn+l + • • • + ttiN Jn- 

For every i, the support Supp (qi) — that is the set of exponents of its non-zero monomials — is 
contained in {A x {0}) U 5 C Z^, where S := {e„+i, . . . , ejv} C and then, by |3^, Prop. 2.12], 

deg^((?i,...,(7,) < Vol((^x {0})U5). 

As we have that 

Vol((^ X {0}) U5) = iV! volRivConv((^ x {0}) U5) ^ n\ volR^Conv(^) = Vol(^) 
(where voIrn and voIr^ denote the standard Euclidean volume forms) we infer that 

5:=(5(/„+i,...,/jv;yo) < Vol(^). 

We conclude that Res^ can be probabilistically computed by means of subroutine Equidim within 
complexity (N n){N dS)'^^^^L{fn+i, . . . Jn) < (n + Vol(^))'^(i). □ 

Remark 4.3 It would be interesting to improve this algorithm in order to compute Res^ without 
any extraneous scalar factor. It would suffice to compute this factor as the coefficient of any extremal 
monomial of Txaj ™s know a priori that the corresponding coefficient in Res^ equals ±1 (^^, 
Ch. 8, Thm. 3.3], see also Cor. 3.1]). 

Example 4.4 We take the following example from Exmpl. 4.. 13]: Set 

A{n, d) := {0, ei, . . . , e„, ei H h e„, 2(ei H h e„), . . . , d (ei, H h e„)} C Z". 

It is easy to check that Vol(^(ri, d)) — nd, and so the previous algorithm computes a sip for ( a scalar 
multiple of) ^esA(n.d) within {nd)'^^-^'> arithmetic operations. 



54 



4.2 Generic over-determined systems 



Our last application concerns the computation of the unique solution of a generic over-determined 
system. 

Let /o, • • • , fn G Q[xo, ■ ■ ■ , Xn] be homogeneous polynomials of degree d. The associated equation 
system is generically inconsistent, where generically means if and only if the vector of the coefficients 
of the polynomials does not lie in the hypersurface V(ReSn^d) C defined by the classical 

resultant Res„_d of n + 1 homogeneous (n + l)-variate polynomials of degree d, and N := (''^") — 1. 
Now assume that the system is consistent. In this case the system is said to be over- determined, in 
the sense that its solution set can be defined — at least locally — with less equations. 
Under this condition the system has generically exactly one sol ution, which is a rational map of 



the coefficients of the polynomials /o, . . . ,/„ (see Corollary 4/7 below). A natural problem is to 
compute this rational parametrization. In what follows we show that this parametrization can be 
easily derived from the resultant, and therefore can be computed with our algorithm. 

In fact we treat the more general case of an over-determined linear system on a variety. The following 
result seems to be classical. However we could not find a proof in the existing literature, so we provide 
one here. 

Lemma 4.5 Let V C P" be an equidimensional variety of dimension r definable over Q. Let 
J-'v{Uq, . . . , Ur) be a Chow form of V , and let u :— (mq, . . . ,Ur) € V{J-v) C (P")''+i be such that 
dJ-y /dUigjg{u) ^ for some < io < r, < jo < n. For < i < r, let Li{Ui,x) := U^xq -\- 
■ • ■ + UinXn denote the generic linear form associated to Ui. Then V fl V{Lq{uq, x), . . . , Lr{ur, x)) 
consists of exactly one element ^{u), and 

Proof.- As the formula stated by the Lemma is invariant by linear changes of variables, we can 
assume w.l.o.g. that no irreducible component of V is contained in any hyperplane {xj = 0}, 
< i < n. 

For < i < r we set ti{x) := Li{ui, x) = UiQXQ + ■ ■ ■ + Ui„Xn G C[a;o, . . . , Xn] for the linear form 
associated to Ui e C"+^. Then V fl V{io, ■ ■ ■ ,(-r) 7^ because of the assumption J-viu) = 0. Let ,f 
be a point in this variety. Suppose 7^ 0. 

Set V'"^^ C A" for the image of V under the rational map '0 : P" A" defined by (xq : • • • : 
Xn) ^ {xi/xq, . . . , Xu/xq). Let T := {Tq, . . . , T^} be a group of r + 1 additional variab les, and let 



P :— PyaS & QlUjlT] be the characteristic polynomial of V^^, as defined in Subsection 3T. Then, 
for < j < n, 

The first equality was shown in Lemma |3.l| , while the second follows directly from formula in 
Subsection |3.l| , and the fact that Li{ui, ^) = for < i < r. 

From Identity (12) and the assumption dJ-y /dUigjg{u) 7^ 0, we infer that dTy /dUigo{u) ^ and 



iu). Therefore 



This shows in particular that £0, ■ ■ ■ ,(r have exactly one common root in \ {xq = 0}. Moreover, 
as the formula for the coordinates of ^ does not depend on the chosen affine chart, we conclude that 
^ is the only common root of £q, . . . , £r in V . □ 
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By the way, we point out a mistake in Ch. 3, Cor. 3.7]. This Corollary would imply that the 
formula of Lemma 4.5 holds in case is a simple common root of £o,. . . ,£r in V. Denoting by 
Ov.^ the local ring of V at ^, this is equivalent to the fact that Ov,£,/{£o, ■ ■ ■ ,£r) = C. 

The following counterexample shows that this is not true: let 

Fit,x,y):=x^{x + t)-ty^ eQ[t,x,y] , C:^V{F)(lP^ 

C is an elliptic curve with a node at (1 : : 0). The linear forms £o '■— io((0 : 1 : 0),{t : x : y)) — x 
and ii := Li{{0 : : 1), {t : x : y)) = y have a single common root (1 : : 0) in C, which is a simple 
root of £o, £i in C. 

On the other hand, as C is a hypersurface, J-c = F{Mq, —Mi, M2), where Mj denotes the maximal 
minor obtained by deleting the (j + 1) column of the matrix 

/ C/00 C/01 U02 \ 
V Uio Un U12 J ' 

A straightforward computation shows that dJ-'c/dUij{{0 : 1 : 0), (0 : : 1)) = for every 

The proof given in [|9) is based on the biduality theorem and on Cayley's trick, and it holds in 
case V is smooth, and in case u — {uq, . . . , Ur) does not lie in the singular locus of the hypersurface 
V{J-v)- This last condition is equivalent to ours. 

Let V C P" be an equidimensional variety of dimension r, and set Vly '■— V{J^v) C (P")''+i for the 
set of (coefficients of) over-determined linear systems over V . As J-y is squarefree and each of its 
irreducible factors depends on every group of variables, 

for < I < r. Then 6^ := fiy \ V{d!Fv /dUiQ, . . . , dTy /dUin) is a dense open set of fiy and so 

( c>T FiT \ 



is a rational map well-defined on 80 U • • • U 



Now let V C (P")g be an arbitrary variety of dimension r, and let V = U ■ • ■ U Vq be its 
equidimensional decomposition. In what follows, for sake of clarity, we keep the same notations as 
previously for different objects sharing analogous properties. 
Set (again) ily C (P")'-+i for the set 

nv - {{Uo, ...,Ur)e (P")''+l -.B^eV / Loiuo, = 0,..., LriUr, = 0} 

of generic over-determined linear systems over V, which is a quasi-projective variety of codimension 
1 in (P")'-+i. 

For every < fc < r, let Qv^ be the set of the coefficients of r -I- 1 linear forms which have a common 
root in Vk ■ If J-^Vi. is a Chow form of Vk , we have that 



0<io<---<ik<r 

and, therefore, has codimension at least 2for0<fc<r — 1. 
Let &i := \ vldTvJdU.o, dTyJdU.n) for i = 0, . . . , r. 

Then every over-determined linear system over V with coefficients in the open set (8o U • ■ • U Qr) \ 
(f^Vb U ■ ■ • Ur2y^_i ) of riy has a unique solution in V which, in fact, lies in Vr- As before, this solution 
can be given by the rational map 



As an immediate consequence of Theorem 1 and Lemma 4.5 we obtain: 
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Corollary 4.6 Let /i, . . . , /s, g G Q[a;o, . . . , Xn] be homogeneous polynomials of degree hounded by 
d encoded by straight-line programs of length bounded by L. 

Set V :— V{fi, . . . , fs) \ V{g) C P" for the guasi-projective variety { /i = 0, . . . , /s = 0, 5 7^ } and 
let V ^ Vq U ■ ■ ■ U Vn be its minimal equidimensional decomposition. Let S := S{fi, . . . , fs] g) be the 
geometric degree of the input polynomial system. 

Then there is a bounded probability algorithm which computes (sip's for) the coordinates of the 
rational map defined above within (expected) complexity s{nd6)'~'^^^L. 

The previous result can be applied directly to compute the solution of a generic over-determined 
system of n + 1 homogeneous polynomials in n + 1 variables of degree d by means of ReSra,^: 

Corollary 4.7 Let u = (uq, . . . , «„) G (P^)«+i where N ('^+") - I, and forO<i<n, set 

fi • — ^ ^ Uia X . 

\a\=d 

Assume that KeSn, d{u) — and that dllesn,d/ dUi„ i3{u) ^ for some < io < n, (3 = (/3o, ■ • • , Pn) G 
(No)"+i with \f3\ ^ d. 

Then V^(/o, ■ • ■ , fn) consists of exactly one element £,{u) G P", and 

I ( 9Res„ d d^c'&n.d , N 
= 7577 W : ■ • • : 7^77 (") 

\Ct^io,(<i-l)ej+eo C^i^ _ (d_ i) e^. +e„ 

for any < j < n such that /3j =/= 0. 

Proof.- From Lemma 4^ applied to the Veronese variety V{n,d) C (see Section [4.1.l| ) we have 
that V{fa, . . . , /„) has only one point ^(u) and that 

«(«r)H- = (^(")) ■ 

V OUioa J \a\=d 

Let /3 = (/3o, • ■ • , (in) be such that \(3\ — d and dReSn,d/ dUigi3{u) ^ 0, and let < < n be such that 
Pj 7^ 0. The previous identity implies that ^ G {xj ^ 0}. Then 

J./ N //-d—1 /- /-d—i \ ( 5Res„ (i 9ReSn d / x 



□ 



As an immediate consequence of this result and Proposition 4.1 we obtain: 

Corollary 4.8 Let notation be as in Corollary [4. '<[ . Then the rational map (P^)*^' ' " F", u 
^{u) can be (deterministically) computed within complexity {nd"')'-'^^K 
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